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We study observables in a conformal field theory which are very closely related to the 
ones used to describe hadronic events at colliders. We focus on the correlation functions 
of the energies deposited on calorimeters placed at a large distance from the collision. 
We consider initial states produced by an operator insertion and we study some general 
properties of the energy correlation functions for conformal field theories. We argue that 
the small angle singularities of energy correlation functions are controlled by the twist of 
non-local light-ray operators with a definite spin. We relate the charge two point function 
to a particular moment of the parton distribution functions appearing in deep inelastic 
scattering. The one point energy correlation functions are characterized by a few numbers. 
For M = 1 superconformal theories the one point function for states created by the R- 
current or the stress tensor are determined by the two parameters a and c characterizing the 
conformal anomaly. Demanding that the measured energies are positive we get bounds on 
a/c. We also give a prescription for computing the energy and charge correlation functions 
in theories that have a gravity dual. The prescription amounts to probing the falling 
string state as it crosses the AdS horizon with gravitational shock waves. In the leading, 
two derivative, gravity approximation the energy is uniformly distributed on the sphere 
at infinity, with no fluctuations. We compute the stringy corrections and we show that 



they lead to small, non-gaussian, fluctuations in the energy distribution. Corrections 
the one point functions or antenna patterns are related to higher derivative corrections 
the bulk. 



1. Introduction 



In this paper we consider conformal field theories and we study physical processes 
that are closely related to the ones studied at particle colliders. In some sense we will 
be studying "conformal collider physics". We consider an external perturbation that is 
localized in space and time near t ~ x ~ 0. This external perturbation couples to some 
operator O of the conformal field theory and produces a localized excitation in the con- 
formal field theory. This excitation then grows in size and propagates outwards. We want 
to study the properties of the state that is produced. For this purpose we consider ide- 
alized "calorimeters" that measure the total flux of energy per unit angle far away from 
the region where the localized perturbation was concentrated. As a particular example 
one could have in mind a real world process e + e~ — > 7* ^hadronsS, where we produce 
hadrons via an intermediate off shell photon. We can treat the process to lowest order in 
the electromagnetic coupling constant and to all orders in the strong coupling constant. 
The QCD computation reduces to studying the state created on the QCD vacuum by the 
electromagnetic current jg m . From the point of view of QCD this current is simply a global 
symmetry. In this case the theory is not conformal, but at high enough energies we can 
approximate the process as a conformal one to the extent that we can ignore the running 
of the coupling and the details of the hadronization process. In this paper we will analyze 
similar processes but in conformal field theories. 




Fig. 1: A localized excitation is produced in a conformal field theory and its 
decay products are measured by calorimeters sitting far away. 



For early work on the applications of scale invariance to strong interactions and, in particular, 
e + e~ collisions, see 
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Our goal is to describe features of the produced state. For example, at weak coupling 
we expect to see a certain number of fairly well defined jets. At strong coupling we expect to 
see a more spherically symmetric distribution 0,^,3 . We need suitably inclusive variables 
which are IR finite. In QCD this is commonly done using inclusive jet observables , see 
for a review. In this paper we study a particularly simple set of inclusive observables 
which are the energy correlation functions, originally introduced in [7[]. They are defined 
as follows. We place calorimeters at angles 9\ , • • • , 6 n and we measure the total energy per 
unit angle deposited at each of these angles. We multiply all these energies together and 
compute the average over all events. These are also inclusive, IR finite observables which 
one could use to study properties of the produced state. Energy correlation functions for 
hadronic final states have been measured experimentally and they are one of the ways of 
making precise determinations of a s (see |§ for example). 

A nice feature of energy correlation functions is that they are defined in terms of 
correlation functions of local gauge invariant operators. They are given in terms of the 
stress tensor operator || . More precisely, consider the expression for the integrated energy 
flux per unit angle at a large sphere of radius r 

/oo 
dtn'T^t.rn 1 ) (1.1) 
-oo 

where n l is a unit vector in R 3 and it specifies the point on the S 2 at infinity where 
we have our "calorimeter". If we integrate this quantity over all angles we get the total 
energy flux which is equal to the energy deposited by the operator insertion. Energy 
correlation functions are defined as the quantum expectation value of a product of energy 
flux operators on the state produced by the localized operator insertion 

(m) ■ ■ ■ Wn)) = jo\ow\o) (L2) 

where O is the operator creating the localized perturbation. Note that the operators are 
ordered as written, they are not time ordered. Notice, also, that the expectation values 



in the left hand side of ( |1.2|) are defined on the particular state created by the operator 
O and they are not vacuum expectation values. The energy operators are very far away 
from each other and they commute with each other. This will become more clear below 
when we think of the operators as acting on null outgoing infinity, sometimes called J~ + . 
Of course, we usually think of the energy deposited at various calorimeters as commuting 
observables, since we measure them simultaneously. Notice that when we compute an n 
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point function we place calorimeters at n points but we also allow energy to go through 
the regions where we have not placed calorimeters. 

In this paper we will assume that we have a conformal field theory. There are several 
motivations for doing so. First, the conformal case is simpler because it has more symmetry 
and, at the same time, it allows us to consider theories that are strongly coupled. There 
are some interesting statements that can be made using conformal symmetry. Second, we 
could have a theory for new physics beyond the Standard Model which is conformal, as in 
the Randall-Sundrum II []10| or the unparticle scenarios, or approximately conformal, 
as in the "hidden valley" scenario [fL2| . One would like to describe the events in these 
theories. In order for energy correlations to be observable to us we need some way to 
transfer the energy from the new sector back to the standard model, as in |TI|. Depending 
on the details, this conformal breaking and conversion process might or might not destroy 
the energy correlations one computes in the conformal theory. We will not discuss this 
problem here. A similar issue arises in QCD. For a sample of references on the influence 
of hadronization on energy correlations for QCD see [|T|,|T3||. The final motivation is a more 



theoretical one, which is to understand better the AdS/CFT correspondence [14jl5 , 16 



Energy correlations are natural observables on the field theory side which one would like 
to understand using gravity and string theory in AdS. We will see that on the gravity side, 
energy correlations translate into the probing of a string state, created by the localized 
perturbation, with a gravitational shock wave as it falls into the AdS horizon. Thus, the 
problem becomes a high energy scattering calculation in the bulk. 

This paper is organized as follows. In section two we make some general remarks 
on energy correlation functions in conformal field theories. By making conformal trans- 
formations we can picture the problem in various ways. We also make some remarks on 
the small angle behavior of the correlators when two of the energy operators come close 
together. We point out that this small angle behavior can be analyzed by means of an 
operator product expansion which involves non-local light-ray operators which are closely 
related to the ones that appear in the discussion of deep inelastic scattering. We also 
relate a moment of the deep inelastic cross section, or parton distribution function, to a 
particular energy two point correlation function. Finally, we consider the general form of 
the energy one point function (£(#)) and relate it to vacuum expectation values of three 
point functions. 

In section three we study conformal field theories that have a gravity (or string theory) 
dual and we describe a prescription for computing the energy correlation functions. The 
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procedure amounts to taking a "snapshot" of the wavefunction of the state produced by 
the operator insertion. In the gravity approximation we find that the energy correlation 
functions are perfectly spherically symmetric as was expected from the very rapid frag- 
mentation that one expects at strong coupling. This phenomenon was originally analyzed 
in deep inelastic processes in @(see also 3). 

In section four we discuss the leading stringy corrections. They amount to small 
fluctuations in the energy distribution of order 1 / y/X. We also consider these corrections 
for charge correlations which have interesting features in the case that the charges are 
carried by flavor symmetries. Finally, we study the regime where two of the angles come 
close together and find that the result is determined by the energy of peculiar non-local 
string states which are dual to the light-ray operators that appeared in the general field 
theory discussion. These operators have a high conformal dimension at strong coupling 
going like A ~ A 1 / 4 . 

In section five we present a summary, conclusions and a discussion of open problems. 

2. Energy correlations in conformal field theories 

In this section we study energy correlation functions in general conformal field theories. 
The discussion in this section is valid for any value of the coupling. 

2.1. Energy correlations in various coordinates systems 

The goal of this subsection is to think about energy correlations in various coordinate 
systems in order to make manifest its various properties and also in order to simplify later 
computations. 

It is interesting to take a step back and think about the energy density as follows. For 
any generator, G, of the conformal group there is an associated conformal killing vector 
(,q (x^ — > + (q). The associated conserved charge can be written as the integral of 
a conserved current, constructed by contracting (q with the stress tensor, over a spatial 
hypersurface 

Qg= [ *±3G , 3^ = T^C G (2.1) 

where the normalization of the stress tensor is chosen so that T^ u = — ^p^pF- This 
expression of the charges is covariant under conformal transformations. It is also invariant 
under Weyl transformations of the four dimensional metric! g^ u — > VL 2 g^ v , T^ v — ► Q~ 2 T flu . 

2 We are ignoring the conformal anomaly since it only contributes as a c-number, independent 
of the quantum state of the field theory. 
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It is convenient to understand clearly the symmetries of the problem. We are interested 
in measuring the flux of energy at large distances. Thus, we focus our attention on the 
boundary of Minkowski space R 1,3 . The conformal generators that leave the boundary 
fixed are the dilatation and the Poincare generators, including the translations and 
the 50(1,3) lorentz transformations. In other words, we have the whole conformal group 
except the special conformal transformations. In order to see that the large r limit in ( |1 . 1| ) 
is well defined, and also to gain some more insight into the problem, it is convenient to 
perform a conformal transformation from the original coordinates to new coordinates 
y M . The new coordinates are such that the future boundary of the original Minkowski 
space is mapped to the null surface y + = 0. The explicit change of coordinates isi 

y + = -^, y"=x--^±^, y 1 ^^, y 2 = ^ (2.2) 

•As •A,' lAy •A,' 

where y ± = y° ± y 3 , and similarly for x ± . The inverse change of coordinates is given by 
the same expressions with x <-> y. The advantage of the new coordinates is that now the 
energy is expressed in terms of an integral over the surface at y + = and we do not have to 



take any limit, such as the large r limit in (1.1). Actually, to be more precise, the surface 



y + = corresponds to the future lightlike boundary of Minkowski space. The energy 
correlation function (|1 . 1|) involves an integral over the past and the future boundaries of 
Minkowski space. However, in the physical situation we are interested in, where we have 
the vacuum in the past, there is no contribution from the past light-like boundary and we 
can focus only on the future boundary. Of course, one could also directly define the energy 
flux operator in terms of an integral over only the future boundary. 

In order to switch between different coordinate systems it is convenient to think about 
R 1 ' 3 as follows. We introduce the six coordinates Z M subject to the identification Z M ~ 
XZ M and the constraint^ 

-(Z" 1 ) 2 - (Z ) 2 + (Z 1 ) 2 + (Z 2 ) 2 + (Z 3 ) 2 + (Z 4 ) 2 = (2.3) 

The usual coordinates on R 1,3 are projective coordinates x^ = z -^+ Z 4. , H = 0,1,2,3. 
The metric induced on this surface, ( |2.3| ), by the i? 2 ' 4 metric is fixed up to an overall 
^-dependent factor. We can choose a metric by choosing a "gauge condition" such as 



3 



This type of coordinates has also been studied in [17]. 



4 Note that Z 1 is the "minus one" component of the vector Z and it does not denote the 
inverse of Z. Hopefully, this notation will not cause confusion. 
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Z 1 + Z 4 = 1. Different "gauge conditions" lead to metrics that differ by a Weyl rescaling. 
The coordinates y^ in (|2.2|) correspond to the choice 



y° z° + z s: y3 z° + z 3 ' yl z° + z 3 ' y2 z° + z 3 ( 2 ' 4 ^ 



In fact, using ( p.4| ) and (|2.2|) we can easily go between the two sets of coordinates. We 
have dx^dx^ = ^ y +^2 ■ We also clearly see that (|2.2|) amounts to a ^ rotation in the [-1,0] 
plane and in the [4,3] plane of R > , which is an element of the conformal group. The 
boundary of Minkowski space is the null surface given by Z _1 + Z 4 = 0. We can think of 
the various generators of the conformal group as the antisymmetric matrices M^- MN ^ which 
generate the transformations 5Z N = M^ NM ^Z^,^. Defining Z ± = Z~ x ± Z 4 , we can see 
that all the generators that leave the surface Z + = invariant are all the ones with no + 
index plus the generator Aft" 1 K In this language the four momentum generators in the 
x coordinates correspond to ilft - ^, /j, = 0, 1, 2,3,4. These generators have a particularly 
simple form at Z + = 

P »~ Z ^- Z -J^ — > P »\z + =o~ Z ^ (2-5) 

(note that Z- = —Z + /2). Since the Killing vectors are all proportional to each other, then 
all four generators involve a single component of the stress energy tensor. Using ( |2.1| ), ( |2.4| ) 
and (^3|) we can write 

P*+Pl= j dy 1 dy 2 £(y lj y 2 ) 

P° ~ p x=j dyidy 2 {y\ + y 2 2 )£(yu 2/2) 

Pi = J dy 1 dy 2 y 1 £(y ll y 2 ) (2.6) 

p x = J dyidy 2 y 2 £{y 1 ,y 2 ) 

/oo 
dy-T__(y-,y+ = 0,y\y 2 ) 
-00 

We see that they are all determined by T__ thanks to the simple form of the generators at 
Z + = (|2.5|) . The conclusion is that we are computing correlation functions of T__ and 
these determine all the components of the energy and the momentum. These expression 



This 50(2,4) manifestly invariant formalism has also been studied recently in [15]. 
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have the advantage that no limit is involved but they have the disadvantage that the 
SO (3) rotation symmetry is not manifest. Since no limit is involved, it is clear that the 
expectation values of ( |2.6| ) will be finite. In fact, we are considering an external operator 
insertion which is localized in x space. This implies, in particular, that it is localized 
near x + ~ so that it is far enough from y + = which is the point where we insert the 
operators ( ^6|) . 

We should note that the dilatation symmetry of the original coordinates x M — > Ax M 
becomes a boost in the y + ,y~ plane in the y coordinates ( |2.2| ). Similarly the dilatation 
transformation in the y variables becomes a boost in the x^ 1 plane. 





(a) 



Fig. 2: (a) Penrose diagram of flat Minkowski space. The doted line is a surface 
at constant r where we measure the energy flux. In the large r limit this becomes 
the light-like boundary, J~ + , of Minkowski space. We consider only the future part 
of the boundary. The semicircle represents a localized operator insertion. In (b) 
we extend the coordinates to the conformal completion of Minkowski space, which 
gives us S 3 x R. The future boundary of the original space is simply the light-cone 
of the point at spatial infinity, i° . 



An alternative point of view is the following. We write the original coordinates as 

" -dt 2 + dr 2 



ds A = -dt* + dr z + r 2 d£l 2 



2.7) 



The original metric and the bracketed metric in ( |2.7| ) differ by a Weyl transformation, 
but such a transformation leaves the physics of the CFT invariant. So we can view our 
CFT as defined on an extremal black hole: AdS 2 x S 2 . Then, the boundary of Minkowski 
space corresponds to the black hole horizon situated at t, r = oo. We see that we can view 
our measurement as one done at the horizon of an extremal black hole. (Of course we 



8 



can also consider other coordinates related by Weyl transformations which would suggest 
other pictures.) By introducing new coordinates we can write the AdS 2 metric in Q2.7| ) as 

l9 —dt 2 + dr 2 —dr 2 + da 2 sinr sincr . „. 

ds 2 = = — , t = , r = (2.8) 

r z sin a cos r + cos a cos r + cos a 

The horizon is at r + = r + a = n. We also define r~ = r — a. We can then write the 
generators ( |2.6|) as 

P° = J dfl 2 £(n) 

P l =J dQ 2 n l £(n) (2.9) 
£{n) =2 [ dr~ ( cos % 



2 , 

where ri 1 is a unit vector in R 3 and specifies a point on S 2 . In these coordinates the SO (3) 
rotation symmetry is manifest. The fact that the energy flux and the momentum flux is 
related to the same operator, T__, is indeed what we would naively expect in a theory 
of massless particles. Namely, if at some point of the sphere we have energy £(9) then 
we have momentum P l = ri l £{9). Here we have shown that this also holds for a general 
interacting CFT. This is due to the simple form of the Killing vector (|2.5| ) at Z + = 0. 

Note also that the SO(l,3) Lorentz symmetry acts on the 2-sphere as the SL(2,C) 
group of conformal transformations of S 2 . Our problem however, does not reduce to 
computing correlators in a 2d CFT, since the state we are considering breaks the SL(2, C) 
invariance. Under these transformations the operator £ transforms as a dimension three 
operator. The easiest way to see this it to recall that these SL(2, C) transformations are 
the ordinary Lorentz transformations of the original coordinates. In particular we have 
seen that x ± boosts become dilatation operators in the y variables. In those variables it is 
clear that J dy~T__ has dimension three. In particular, one can find the relation between 
the operator £(y 1 ,y 2 ) which is defined on a plane to the one on the sphere, £(n), by 
following the coordinate transformation between the plane and the sphere at Z + = y + = 

sinfle^ 6 itD 

yi + iy2 =77- = tan -e v 

(l + cos0) 2 



(1 + COS ft) 2 

£(yi, y 2 ) =n- 3 £(n) = (1 + cos#) 3 £(n) 



Physically, we expect that our idealized calorimeters will measure positive energies. 
Therefore, the expectation values of £ (n) should be non-negative. In quantum field the- 
ory the expectation value of the stress tensor can be negative in some spacetime region. 
However, in our case we are integrating the stress tensor along a light like direction. In a 
free field theory one can show that the expectation value 



/ 



dy-(T__)>0 (2.11) 



is positive on any state |19|]ll. We expect that the same should be true in an interacting 
field theory. In appendix A we recall the argument in free field theories and give a hand- 
waving argument suggesting that this should be true in general. We will later see that this 
condition implies interesting constraints on certain field theory quantities, so it would be 
nice to be able to give a more solid argument for the positivity of ( j2.11p than the one we 
give in the appendix. 

Notice that the energy flux operators S(9) commute with each other since operators 
at different values of 9 are separated by spacelike distances. This is most clear when we 
express the operators in terms of the y coordinates as in ( ^.6| ). Thus, we can certainly 
consider the probability that we measure specific energy functions S(6) = f(6) and derive 
the probability functional that governs the process. Once can also impose some cuts on the 
energy distribution and compute such probabilities. This is done when jet cross sections 
are computed, as in ||. In fact, a specific Feynman diagram with n particles coming out 
at angles 6*i, • • • , 6 n gives a contribution to the case where the energy function f{6) is a 
delta function localized at these points. The energy correlation functions we have defined 
correspond to average energies where we also allow extra particles that come out and do 
not go into the calorimeters we are choosing to focus on. 

Besides putting a detector at infinity that measures energy we can also put a detector 
that measures charge. In that case we have the charge flux operator 

/oo 
dtn l j t (t,rn) (2.12) 
-oo 

where j is the current associated to a global U(l) symmetry of the field theory. In the 
coordinates (|2.2|) this becomes Qiy 1 ^ 2 ) = j dy~ j-(y~ ,y + = O,?/ 1 ,?/ 2 )). Under SL(2,C) 



6 The curved space analog of this condition has also been explored for free fields in curved 
space, since it plays a role in proving singularity theorems in general relativity. 
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transformations Q transforms as a field of conformal dimension two. We can similarly 
compute energy and charge correlation functions. One can also easily consider non-abelian 
global symmetries, and measure the components of various charges, as long as we do not 
put two charge insertions at the same point. 

Now let us make some remarks on the operator ordering. Since the energy flux op- 
erators commute with each other for different 9, then, it does not matter how they are 
ordered. However, it is important that they are inserted between the operator, O, that 
creates the state and the one annihilating it, as in ( |1.2| ). This is the standard ordering 
when we compute expectation values. If we use perturbation theory to compute them it 
is important that we do not use Feynman propagators since those are for time ordered 
situations. However, to do perturbation theory it is very convenient to use Fenymann 
propagators. In such a case we have to be careful to remember that we should use the in- 
ZD| , |2I|] formalism to evaluate the expectation value. This consists in choosing a contour 



m 



that starts with the initial state, goes forward in time to the times where the stress tensor 
operators are evaluated and then goes backwards in time. 

In a conformal field theory we could also consider the following. Minkwoski space can 
be mapped to a finite region of R x S 3 . In fact, Rx S 3 can be split into an infinite number 
of regions, each of which is mapped to Minkowski space. In that case we can consider one 
of the regions as the original Minkowski space and the region immediately to the future 
as the region parametrizing the part of the Schwinger-Keldysh contour that goes back in 
time, as long as we transform the wavefunction of the in state in the bra appropriately. 
We found this picture useful for gaining intuition, but not particularly useful for doing 
computations. 

2.2. Small angle singularities and the operator product expansion 




(a) < b > 



Fig. 3: (a) Singularities in the energy correlation functions arise when we place 
two calorimeters very close to each other, at a small angle 6. (b) At the level of 
Feynman diagrams such singularities come from colinear radiation. 
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The energy correlation functions develop singularities when two of the energy opera- 
tors are evaluated at very similar angles 9\ ~ #2, see fig. 3a. Such singularities are related 
to collinear radiation. To leading order in the gauge theory coupling A = 47ra s N the 
leading singularity goes like £(#i)£(#2) ~ an d h comes from a Feynman diagram 
like the one shown in fig. 3b. 

It is clear that such a limit should be characterized by some sort of operator product 
expansion. In this section we will make some remarks on the type of operators that appear 
in this expansion. 

It is simpler to think about the problem in the y M coordinates introduced in ( |2.2j ). We 
should, then, compute the OPE of operators of the form 

£(y\y 2 )£(0,0)~ J dy-T__(y-,y+ = 0,y) J dy'~T — (y'~ ,y + = 0, 0) (2.13) 

The two operators are sitting at two different points in the transverse directions. We have 
set one at zero for convenience and the other at y = (2/1,2/2)- Note that the distance 
between the two stress tensor insertions is \y\ irrespective of the values of y~,y' ■ This 
distance is spacelike, so one expects to be able to perform an operator product expansion 
when ?/* — > 0. Nevertheless, since the two stress tensors are sitting at two very different 
points in the y~ directions, the operators appearing in the OPE are not local operators. 
To leading order the operator is specified by two points that are light-like separated |22T| . 
Such operators are useful for thinking about many high energy processes in QCD [p3| , [2~^ , p5| . 
They are sometimes called "string operators" or "light ray" operators. Various "parton 
distribution" functions are defined in terms of matrix elements of such operators, see |24| 



for example. It is important to note that these operators are non-local along one light-like 
direction but they are perfectly local in all remaining three directions. 

In order to characterize these non-local operators it is useful to label them according 
to their transformation properties under the conformal group |26| (see [|27j for a review). 
Let us define the twist generator to be T = A — j, where j is the spin (really a boost 
generator) in the j/ + , y~ plane@. More explicitly, the twist transformation is (r/ + , y~ , y) — ► 
(A 2 y + , y~, Xy). The spin is the transformation (y + ,y~,y) — > (rjy + ,rj~ 1 y~ ,y). As it is 
well known, at zero coupling, one can consider twist two operators which correspond to 



7 Note that we define j to be the spin in the y + ,y plane only, not the total spin. The spin in 
the transverse directions is another generator which does not appear in the definition of the twist. 
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primary operators of higher spin. For example, if we have a scalar field, <f>, in the adjoint 
representation, then we can schematically define the operators 



Uj = Tr[<ffiL<f>] (2.14) 

This is schematic because there is a precise combination of derivatives that makes it a 
conformal primaryi. Such conformal primaries exist only if j is even. One is sometimes 
interested in extending the definition of such operators to generic, real or complex, values 
of j. This problem was considered in detail in ||26|| . There, it was found that one could 
start with the operators 

U(y-,y'-) = Tr[<j ) (y-)W(y-,y'-) ( p(y'-)} = Tr[<p(y-)Pefy- A <P(y'~)] (2.15) 

where W is an adjoint Wilson line along a null direction. All operators are inserted 
at the same values of y + ,y x and y 2 (but of course, at different values of y~). We can 
also replace by a fermion or a gluon operator F-i. Under twist transformations y~ 
remains invariant but the transverse coordinates are rescaled. In the quantum theory 
this scaling transformation mixes the operator ( |2.15|) with operators with other values of 
y~,y' ■ By thinking about the action of the collinear conformal group (the SL(2,R) set 
of transformations of x~) it is possible to diagonalize the action of the twist generator. To 
leading order the operators are diagonalized by considering suitable combinations of these 



light-ray operators |2q , p7 | . These operators are labeled by their center of mass momentum 
k- along the y~ direction and their spin. For our purposes we will be interested only in 
operators which are integrated over the center of mass position along the y~ direction so 
that they carry zero momentum along y~ . In that case the operators of arbitrary spin 
constructed from scalar fields can be written as 

poo 

u 



<J-1 



/ dy~ / -— Tr[cf>(y- + u)W{y~ +u,y~- u)<j>{y- - u)] (2.16) 

J-oo JO u 



The subindex of U denotes the total spin and j denotes the spin before we do the y~ 
integration. This is an expression that makes sense for arbitrary complex values of j. When 
j approaches an even integer we find a pole in j coming from a logarithmic divergence in 
the integral at small u of the form J — . The coefficient of this divergent term contains the 



The precise form is Uj = J2i =0 rfc. ( r-M n-2 Tr[(d k (f))d J r k (j)}, where <f> is a scalar field @|9],|3(|. 
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ordinary local operator ( 2.14 ), see [26,2?| for more details. There are similar expressions 



for operators constructed from two fermions or two Yang-Mills field strengths. One can 
compute the value of the twist for these operators and one finds [^6[ r(j) = 2+j(j), where 



7(j) is the anomalous dimension. One can also consider higher twist operators which 
contain more field insertions or extra derivatives with respect to the transverse direction 
or y + . In that case, in order to diagonalize the matrix of anomalous dimensions, it is not 
enough to give the total spin of the operator. Nevertheless, this can be done, see |]2^ . 
The OPE has the schematic form 

£(y)£(0)~ I dy-T__(y-,y) j dy'-T^y 1 ' , 0) ~ £ \y\ T "~ 4 ^-1,^3 (2.17) 

n 

where the sum is over all operators which are local in y + ,y, but not necessarily local in 
y~ , which have total spin j — 1 = 2, (or j = 3) and twist r n . The spin is determined since 
the total spin of the left hand side is one for each of the two energy insertions. Equation 
( |2.17|) is schematic because we have not explicitly indicated the fact that the operators in 
the right hand side could carry spin in the transverse directions. A more precise expression 
has the form 

£(y)£(0)~~£y^--V fc) iyr^^ (2-18) 

k,n 

where we have now considered operators that carry spin in the transverse directions, the 
indices ii, • --ifc are symmetric and traceless. 

Among the operators which have twist two at zeroth order there are only a few that 
have j = 3. For example, in QCD there are only two, a bilinear in fermions and a bilinear 
in the gluon field strength. Thus, for the given spin we are considering (j = 3) we will 
have to diagonalize a finite matrix of anomalous dimensions. 

In summary: The small angle behavior of the energy correlation functions is deter- 
mined by the spin j = 3 non-local operators that appear in the OPE 

(£(e 1 )£(6 2 ) \0i2\ T "- 4 (^3-iA^) ■ ■ ■) (2.19) 

n 

where the dots denote other energy insertions and \9\ is the angle between the two energy 
insertions that are getting close to each other. The sum over n runs over all the higher twist 
operators that can appear. We will see that in Af = 4 super Yang Mills these operators 
develop large anomalous dimensions at strong coupling. 
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Note that the spin symmetry in the y + , y~ plane, that we used to select the operators 
that contribute, is the dilatation symmetry of the original Minkowski space. This symmetry 
ensures that the energy correlation functions scale properly as we rescale the total energy 
(or rescale the variables x^). In other words, there can be no anomalous dimensions under 
total energy rescalings since that would conflict with energy conservation. This is the 
physical reason why we are forced to select particular operators in this OPE. 

In the case of QCD the small angle behavior of energy correlation functions was 



computed a long time ago in [|31],|32| using a slightly different language. They also needed 
to include the effects of the running coupling. 

Let us now turn to the case of M = 4 super Yang Mills at weak coupling. The weak 



coupling computation of the leading twist anomalous dimensions was done in |H3|,[M| (see 



also [ |3"5"| , |3"S| ] ) . We should consider operators which are invariant under all the symmetries 
that leave the particular component of the stress tensor in (|2.18|) invariant. These include 
the 50(6) R-symmetry and a parity symmetry. We can classify the operators according 
to their transformation properties under the SO(2) group that transforms the transverse 
coordinates. All operators are made out of a pair of scalars, fermions or gauge field 
strengths. The local operators with zero transverse spin in SO(2) and spin j (j even) in 



the H — directions are 33 



Tr[(0L(p] , Tr[F_*d*- , Tr[?/T_ cf! V] (2-20) 

Supersymmetry relates these three towers of operators. Since supersymmetry carries spin, 
the various members of the supermultiplet have different spin. However, the anomalous 
dimension for all the members of the supermultiplet is the same and it is given by a function 
which has the weak coupling expansion [j3~3|,P 



iti) = Am-v-m} + --- (2.2i) 



where ip = V (z)/T(z). This was computed also to two and three loops in pi ]. The fact 
that j(j = 2) =0 corresponds to the fact that the stress tensor is not renormalized. 

Since we are interested in operators with a definite spin, we conclude that the three op- 
erators that diagonalize the anomalous dimension matrix are in three different multiplets. 



For spin three operators we have the anomalous dimensions [33| 



n - 2 = 7 (,7 = 3) , r 2 - 2 = 7 (j = 5) , r 3 - 2 = 7 (j = 7) 

t 9 - A t 9 - 1U t 9 - 137A (2 ' 22) 

Tl " 2 "2^' T2 " 2 -12^' T3 " 2 -l20^ 
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where we just gave the first order expression. We see from (|2.22 ) and ( |2.21|) that all three 



anomalous dimensions in ( |2.22| ) are positive and T\ — 2 is the smallest one which will give 
us the leading order singularity. However, for weak coupling all three contributions are 
similar. 

In addition to the operators we discussed, we can also have operators which have 
non-zero transverse spin. At twist two, the only one consistent with the symmetries is the 
spin two operator 

U (u);j =Tr[F_ (l ^r 2 F_ l) } (2.23) 
where the indices i, I = 1,2 are symmetrized and traceless. In the M = 4 theory these 



operators are in the same supermultiplet as the ones considered above [gS[ . For this reason 
their anomalous dimension is also given in terms of the same formula 

fj-2 = 7 (j + 2) , f 3 = 2 + (2-24) 

Thus we expect to have a small angle singularity of the form 

3 

(£(y)E(0) ■ ■ •} ~ Yl \y\~ 2+(Ta ~ 2) Ca(U a ■■■) + y^y l) \y\- 4+(f3 - 2) 5(^i) ■ ■ ■> (2-25) 

a=l 

where the operators U a are the linear combinations that diagonalize the anomalous dimen- 
sion matrix for the operators with zero transverse spin. c a and c are coefficients that can 
be obtained by performing the operator product expansion explicitly. These constants are 
independent of the state for which we compute the energy correlation. Of course, the terms 
(U a ■ ■ •) and (Uun • • •) do depend on the state on which we compute the energy correlation 
function. The coefficients c a , c start at order A at weak coupling since it is easy to check 
explicitly that at tree level there is no contribution to the operator product expansion of 
two energy flux operators. 

In QCD one can do a similar analysis, including the effects of the beta function, see 
32| . In that case, the operator made out from scalars in ( |2.20| ) does not contribute. 

Having done one OPE, we could also do a further OPE of the resulting operator with 
a third energy flux operator. That would give an operator of total spin j — 1 = 3, or 
j = 4, and so on. More generally, we can consider the case where n energy operators come 
close together. If we keep the ratios of angles between these n points fixed, then the small 
angle behavior is given by the anomalous dimension of the operator of spin j = n + 1. 
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The structure of a jet at weak coupling is largely controlled by these operator product 
expansions. 

Note that, at weak coupling, after we consider the effects of the anomalous dimensions, 
the energy correlators have small angle singularities that are integrable. In fact, if we do 
the integral over a small angle 9q of the energy two point function we find schematically 

d 2 es{e)s (o) ~ J B ° d 2 0^x ~ (0 O ) 7 * A (2.26) 

where the anomalous twist of the spin three operator is r — 2 = 7* A and 7* is a numerical 
constant. This expression is schematic because at weak coupling we have to include all 
the terms in (|2.25|). If #0 is fixed and A —>■ then we see that the integral gives a finite 



answer. This is to be expected since the total integral of one energy insertion over the 
whole sphere should give the total energy, independently of A I. The fact that we get a 
finite contribution from this region is consistent with the idea that the energy is going out 
in localized jets. We can also estimate the angular size of jets, by finding a #0 i n such a 
way that we get a fixed fraction, /, of the total energy in the jet. This gives an estimate 
#0 ~ e _c / A , where c depends on /. This was originally discussed in ||, see also |32| for a 
more detailed discussion. 

We could also do an OPE of two charge operators, each of which has spin zero, after 
we integrate the spin one current over y~ ( |2.12|) . In this case we get operators with total 



spin j — 1 = 0, or j = 1. Some of these have negative anomalous dimensions. In fact, we 
expect that charge correlators would be more singular at small angles due to the fact that 
a gluon can create a pair of oppositely charged particles fairly easily and there is no reason 
that we couldn't get a divergence when we integrate the charge correlator at small 9. 

Finally, we should mention that in QCD the energy-energy correlation two point func- 
tion was computed for all angles in0 |7|, f40| , |3^j42l , f43j . It was also compared to experiment 



in ||, where it was used as a way to measure a t 



9 Here we are also assuming that the energies are locally positive. For charge correlators we 
cannot make the same argument because the charge can be positive or negative. 
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The results presented in the following references show some disagreements. For a detailed 



comparison between these results see [39|. We thank S. Catani for pointing this out to us. 
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2. 3. Energy flux one point functions 

In this section we will make some simple and general remarks about the energy flux 
one point function 



(0\ OlS(n)O q \0 ) 

\olo q \o) 



These one point functions are determined up to a few coefficients by Lorentz symmetry, 
even in non-conformal theories. Here we will consider these in the CFT context in order 
to make contact with other results in conformal field theories. 

The energy flux one point function ( |2.27| ) amounts to computing a three point function 
in the CFT. Three point functions in a generic CFT are determined up to a few numbers 
by conformal symmetry |¥|,[I5],^6[ . 

Let us start with the case that we create the external state with a scalar operator 
with energy q and zero momentum. Strictly speaking such an operator is not a localized 
insertion. Thus, more precisely, we will be considering operators of the form 

„2 

a' 



O q = / d^xO{x)e- iqxo exp{- X ° + Xl + / 2 + ^ } , qa > 1 (2.28) 



where the last inequality ensures that the operator is localized, has finite norm and has 
four momentum approximately (p = (q, 0) +o(l/a). In particular we have q° ~ q. Once we 
know this precise form of the operator we see that we can also write it in other coordinate 
systems by performing the suitable conformal transformation and taking into account the 
conformal transformation properties of O(x). 

In what follows we will consider field theory states produced by scalar operators, 
O ~ 5, conserved currents O ~ e^, and the stress tensor, O ~ e ijTij- In all cases we 
consider states with essentially zero spatial momentum as in fl2~28l) . The case where o M is 
a generic four vector can be obtained by performing a simple boost of the configurations 
we discuss. 

In the case that we insert a scalar operator it is clear by 0(3) symmetry that the 
energy one point function is constant on the two sphere. In addition the integral over the 
angles should give the total energy. Thus, for a scalar operator we have 

(S(n)) = £ (2.29) 

Even though we know the answer already, it is possible to do the calculation explicitly 
by writing down the unique general expression for the three point function of two scalars 
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and the stress tensor [HSJ. Its normalization is fixed by a Ward identity in terms of the 
two point function of the two scalars. This Ward identity is another version of the energy 
conservation argument that we used above. Writing down the three point function and 
doing the integrals in the limit we indeed obtain (|2.29| ) . One has to be careful about 
the operator ordering. In appendix C we do this explicitly. 

We now turn to the case where the external perturbation couples to a conserved 
current in the CFT. In that case the operator is given by O ejQ ~ ^j^(q) where & is a 
constant polarization vector. Due to the current conservation condition we can identify 
t v 4- \q^. So we can choose e to point in the spacelike directions. In this case 

0(3) symmetry and the energy conservation condition constrain the form of the one point 
function to 



(0|(e*-jt)g(fi) (j-e)\0) q_ 
[C[U>) (0\(e* ■ ji) (j ■ e)\0) 4tt 



1 + a 2 



\e-n\ 2 1 



W 3 



_Q_ 



l + a 2 (cos 2 9- -) 



[2.30) 



where 6 is the angle between between the point on the 5" 2 , labeled by n l , and the direction 
of the polarization vector e\ 

The fact that we have one free parameter is in agreement with the general analysis of 
the three point function of two conserved currents and the stress tensor. In fact in j|5|| it 



was shown that the three point function is determined by conformal symmetry up to two 
parameters and one of them is fixed by the Ward identity of the stress tensor. 

Note that a 2 in ( |2.30|) obeys a constraint that comes from demanding that the ex- 
pectation value of the energy 8(6) is positive, see (|2.11|) . This condition leads to the 
constraint 

3 > a 2 > -- (2.31) 

This one point function was computed for the electromagnetic current in QCD in 0. 
To first order in a s the result is 



3 9a s 

«2 = -^ + ^ + --- (2.32) 

To the order written in ( |2.32| ) we can approximate the QCD computation by a conformal 
field theory with the value of the coupling set by the energy of the process a s = a s (\q°\). 

In the application to e + e~ collisions that produce a gauge boson which in turns couples 
to a current the polarization vector of the current depends on the polarization states of 
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the e + and e~ as well as the type of gauge boson we are considering (7 or Z, Z' , etc). 
In the case that we consider unpolarized electrons we can express the answer in terms of 
the angle with respect to the beam axis, 9b- (cos 9b = n.z where z is the beam axis). The 
polarization vectors for the current are orthogonal to the beam direction and we should 
average over them. After doing this average, we find that (|2.30| ) becomes 



(fr]) = E 5 (0|(e:-,?t)g(n) (j-e s )\0) q_ 

^ {U)) E s (0|(6*-jt)(i. es )|0) - 47T 



1 + a 2 (^sin 2 9 b - - 



(2.33) 



where we sum over polarization vectors transverse to the beam. For a current that couples 
to free fermions we find the familiar (1 + cos 2 9b) distribution, as we can check from the 
leading order QCD result ( ggg ). 

For a current that couples to free complex bosons of charges q\ and Weyl fermions of 
charges q™^ we get 



a f 2 ree = 3 (2.34) 



where we sum over both left and right Weyl fermions. Note that the case where we only 
have bosons saturates the upper bound in (|2.31|) and free fermions saturate the lower bound 
in (2.31). In fact, going back to ( |2.33| ) we see that we get the well known distributions 



proportional to sin 2 9b or (1 + cos 2 9b) for free bosons and fermions respectively. 

We can consider a similar problem now in an Af = 1 superconformal theory. If the 
current is a global symmetry that commutes with supersymmetry (a non-R symmetry) 
then one can see that a 2 = 0. In a free supersymmetric theory we see from (|2.34|) that the 
bosons and Weyl fermions cancel each other. For an interacting theory this follows from 
the fact that such a current is in the same multiplet as a scalar operator, and for a scalar 
operator we do not have any arbitrary parameters P7 ]. Thus the value of a 2 is fixed by 



superconformal symmetry. However, since we got a 2 = in a free theory, we have a 2 = 
for any global symmetry of a SCFT. 

On the other hand we can get a non-zero value of a 2 if we consider the R currentJHI. 
The R current is in a different supermultiplet. In fact, it is in a supermultiplet with 
the stress tensor. All three point functions among elements of this supermultiplet are 
determined by two numbers, c and a ||47[|. These numbers also characterize the anomalies 
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We thank Scott Thomas for pointing this out. 
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of the R current, which are encoded in parts of the jjj and jTT three point functions 
H,|7|. They also contribute to the conformal anomaly on a general background, 



Tff = ^W^saW^ - j^E , E = R^ Sp R^ Sp - AR^R^ + R 2 (2.35) 

where W is the Weyl tensor and E the Euler density. The c coefficient is the only constant 
that appears in the two point functions of the currents and the stress tensor |[47|| . Thus 
c appears in the part of the three point function that is fixed by the Ward identity. The 
coefficient 0,2 is given by a linear combination of a and c. The particular linear combination 
is independent of the theory. It is fixed by supersymmetry. We can compute the precise 
combination by considering the particular case of free field theories. As we explain in more 
detail below we find that 

(8(6)) = l + 3^^(cos 2 6>- h (2.36) 

C o 

This formula was obtained as follows. We used that for a free supersymmetric theory with 
ny vector multiplets and ns chiral multiplets we have 48a = 9ny + ns, and 24c = 3ny + ns 
4£|49,5(|. The vector multiplet has one Weyl fermion of charge 1 and the scalar multiplet 



in a free theory has a Weyl fermion of charge —1/3 and a boson of charge 2/3. Then 
using ( |2.34j) we obtain ( |2.36| ). Note that even though we used free field theories to fix the 
numerical coefficients, the final result ( |2.36| ) is true for a general interacting M = 1 SCFT. 

In M = 4 super Yang Mills a = c and the result for the one point function is spherically 
symmetric. Of course this is not a surprise since U(l) subgroups of the SO(6) symmetry 
group can also be viewed as global symmetries from the point of view of M = 4 written 
as an M = 1 theory. Thus, in M = 4 super Yang mills the result is independent of the 
coupling. 

The positivity constraint ( |2.31|) , together with ( |2.36| ) gives ^ > a > 0. 

We can also consider the energy one point function in the case that the state is created 
by the stress tensor. As above, we take the momentum of the inserted operator in the time 
like direction. Then the operator that we are considering is characterized by a symmetric 
polarization tensor e u which we take to have indices in the purely spacelike directions by 
using the conservation equations. Since the stress energy tensor is traceless, we also take 
e" = 0. By 0(3) invariance we see that the most general form of the three point function 
is 



(0\e* J T lJ S(9)e lk T lk \0) _ q ° 
(Ole^ejfcTjfclO) ~ 4n 



(2.37) 
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We see that we have two undetermined coefficients. This agrees with the generai analysis 



of stress tensor three point functions in ||45 |, where they found that conformal symmetry 
determines the three point function of the stress tensor up to three coefficients, one of 
which is fixed by a Ward identity. We have chosen the constants in the last two terms in 
( |2.37|) in such a way that the corresponding terms integrate to zero on the sphere. 
By demanding that {£) is positive we get the constraints 

, to Ha . 

(i-4 — - ) > o 

V 3 15 ; ~ 

2(l-|-^)+t 2 >0 (2.38) 

3 , to 2t\ . 

-(1 - — -) +t 2 + t 4 > 

2 V 3 15 ; 

We obtain these constraints by using 0(3) invariance in ( |2.37| ) to set n = z. We then view 
the resulting equation as a bilinear form on the space of e's. This space can be divided into 
three orthogonal parts according to their 0(2) invariance properties (the spin of e along 
the z axis). We have an SO(2) scalar, a vector and a symmetric traceless tensor. On each 
of these subspaces we get each of the constraints ( |2.38| ). Each of this limits in saturated 
in a free theory with no vectors, no fermions or no bosons respectively. The fact that the 
first equation is saturated in a theory without vectors is clear. In that case, if we consider 
a stress tensor insertion with spin +2 in the z direction we cannot have emission of bosons 
or fermions in the z direction due to the orbital angular momentum wavefunctions. It is 
also possible to write a general bound on the two coefficients that appear in the conformal 
anomaly ( |2.35| ), see appendix C. 

In an M = 1 supersymmetric theory we find that 

t 2 = 6(c-a)/c , t 4 = (2.39) 



By requiring that (|2.37|) is positive for all choices of traceless tij we find 



\c > a > °- (2.40) 

Of course c > due to the positivity of the two point functions. The bounds are saturated 
by free theories with only vector supermultiplets (upper bound) or only chiral supermul- 
tiplets (lower bound). It is interesting that the lower bound that we obtain in this way 
is precisely the same as the bound obtained in |51[](see also p2| ) based on causality for a 
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gravity theory that contains only the Einstein term and a R 2 term@. In the theory con- 
sidered in |51]] one would also have £4 = 0, though it is not clear whether it corresponds to 
any dual quantum field theory. Here we have only used general field theory considerations. 

For a non-supersymmetric theory it is also possible to derive a bound from ( |2.38j ) . As 
explained in appendix C we find 

^> a ->\ (2.41) 
18 c 3 

where the lower bound is saturated by a free theory with only scalar bosons and the 
upper bound by a free theory with only vectors. Note that the bound in supersymmetric 
theories ( |2.40|) is more stringent than in non-supersymmetric theories (|2.41| ). Let us also 
add that the results from appendix C also allow us to calculate this bound for Af = 2 
supersymmetric theories. In this case we can obtain the bound by taking the operator O 
to be one of the SU(2) R-symmetry generators and demanding that the energy one point 
function is positive. The result is in this case 

5 - > a - > 1 (2.42) 
4 - c ~ 2 y ' 

This is a smaller window than for the M = 1 expected. The upper bound 

corresponds to a free theory with vector supermultiplets only while the lower bound cor- 



responds to a free theory with hypermultiplets only. This agrees with results in [54 



We can make similar remarks for operators that involve charge correlations. For 
example, we could consider a theory with an SU(2) global symmetry and then select one 
U(l) C SU{2) to form the charge flow operator Q that we measure at infinity. We could 
consider a charged state state created by the current e- J + , where the plus indicates that it 
carries charge plus one. As in the energy correlations the charge correlations have a form 



= (0|e»- J-Q(n)e- J+|0) = J_ 
(0|e* • J- e ■ J+|0) 4tt 



1 + a 2 



n\ 2 1 



W 3 



(2.43) 



Again, the coefficient S2 is related to the fact that there are two possible (parity preserving) 
structures for the three point function of three currents |44]j45[[ . One of them is fixed by 
the Ward identities in terms of the two point functions, a fact we used in (|2.43|) . We note 
that in a supersymmetric theory where these currents are global symmetries 0,2 = 0. One 



can show this as follows. First note from [47] that there is only one parity preserving 



12 In order to see this one has to set \gb — > 9/100 (the bound in |52| ) into the expressions for 
a and c [|| (see eqn. (5.1) of |j§). 



23 



structure for current three point functions in supersymmetric theories. This means that 
the value of 0,2 is fixed by supersymmetry. One can show that is vanishes by computing it 
in a particular theory, such as a free theory or M = 4 super Yang Mills at strong coupling. 

There are some cases where there are parity odd structures that can contribute. Such 
parity odd parts of three point functions are related to anomalies. For example, in the case 
of three currents these are related to the usual anomaly |[44| . Consider the case that we 
have an external state produced by a current and we measure a charge one point function 
distribution. For example, we can consider a U(l) current that has a cubic anomaly. A 
concrete example is the R current in superconformal theories. We consider a state obtained 
by acting with this current on the vacuum and we measure the charge flux far away. We 
find 

?Q(n)e-j\ 0) 
\?-??-j\0) 

Note that this is non-zero only if e is complex. This happens, for example, when we consider 
a circularly polarized state. Then x is the angle between the direction of the spin of the 
current and the calorimeter. This leads to a charge flow asymmetry. Such asymmetries are 
extensively studied in e + e~ collisions that produce a Z boson that then decays. Here we 
are pointing out that the charge flow asymmetries are related to the anomaly. Of course 
the full electroweak symmetry is not anomalous. But if one focuses only in decays of the 
Z into leptons, then the fact that the purely leptonic theory is anomalous leads to the 
charge flow asymmetry. The fact that tree level processes plus unitarity fix the anomaly 
was pointed out in |)5] . 



(Q(n)) = V t A„ =iae jlk e*em k ~ acos X (2.44) 



The three point function of two stress tensors and a current has a term that reflects 
the mixed gravitational anomaly J56|. Consider the case where the stress tensor creates 
the state and we measure the charge distribution of the current that has a mixed anomaly. 
The charge distribution has the structure 

{m} = < '^>^°> = g sagf^ (2 , 5) 

where f3 is related to the anomaly coefficient. For a supersymmetric theory, and for Q 



given by the R-current, we have that j3 ~ (a — c) p8|j5q| . Notice that there is, in principle, 
another tensor structure consistent with 0(3) symmetry that could have contributed to 
fl2.45|) , namely ey^e^e^-n • This term is, however, absent from the three point function 
once conformal symmetry and the Ward identities are imposed [[56] . Thus in a theory with 
a gravitational mixed anomaly there is charge asymmetry for the states produced by the 
graviton. 
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2.4- Relation to deep inelastic scattering 

In this section we explore the relation between the energy correlation functions and 
the deep inelastic scattering cross sections. 

The deep inelastic cross section for the scattering of an electron from a proton can be 
factorized into the electromagnetic process and the strong interactions process. At lowest 
order in the electromagnetic coupling, but exactly in a s , the strong interactions part of 
the cross section can be written in terms of the expectation value of two currents in the 
state of the target (which is traditionally a proton, but can be generalized to any other 
particle) 

= J d 4 ye iqy (p\JV(0)J"(y)\p) = 

=A(*, q 2 /p 2 ) (g»„ - Q -f-)+ ^ 2 (z, f/p 2 ) (p„ + |) (pu + |) (2.46) 

q 2 

where x = , q 2 > 

zp.q 

We are imagining that we have a plane wave state in the y coordinates with timelike 
momentum p 2 < 0. The tensor ( |2.46| ) is nonvanishing only if we create a timelike state 
s = — (q + p) > with the current. For these values of p and q our definition of W^, 
( |2.46| ), coincides with the ordinary one ||57|| , which involves a commutator of the currents. 

We would like to relate these formulas to the ones appearing in the energy correlators. 
Let us consider the charge operator Q evaluated in the y coordinates 



Q(y) = J dy~j-(y + = 0,y~,y) (2.47) 

where y denotes two transverse dimensions. These two transverse dimensions are related 
to the angles on the two sphere by (|2.10| ). In the y coordinates, we can Fourier transform 



this operator. We, then, have something similar to the current appearing above, except 
that the current in ( |2.46|) is in momentum space also in the y + direction. Note also that 
q- = 0, due to the y~ integral in the definition of the charge flux operator Q. Since q- is 
zero, we find that q 2 = (q) 2 > and independent of q + . However x depends on q + since 

1 -2p.q -(p + q) 2 + p 2 + q 2 4p_q + 

- - — 5— 5 ~ — 2 — 1 — y ZA °) 

x q z q z q z 

where the dots indicate terms that are independent of q + . In order to produce the S(y + ) 
that is present in the charge flux operator we need to integrate over q+. This integral 
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translates into an integral over x. The range of integration can be determined by the 
condition that — (p + q) 2 > 0. Thus we end up integrating between x = and x max with 
l/x max = 1 + p 2 /q 2 . In the limit p 2 /q 2 — > we get the usual boundary x = 1. 
We then have the following relation between the two quantities 



J d 2 ye^ p (Q(0)Q(y)) p = ^W—(q+,q- = 0,$p) 

/ — F 2 (x,g /p ) 
Jo x 



(2.49) 



47T 



This is a particular moment of the parton distribution functions. More precisely it is the 

(2) 

moment M| ' . As it is well known, the even moments M^k can be expressed in terms of 
the expectation values of local operators with spins j = 2k ||57|| , via a dispersion relation 

(2) 

argument. In fact, the moments M can also be expressed in terms of the expectation 



values of the non-local light-ray operators with spin j for any j, see |58j for a general 
discussion. 

In ( |2.49| ) the charge correlation is evaluated on a state with definite momentum in 
the y coordinates. This implies, in particular, that the charge two point function is also 
translation invariant in the transverse space, p (Q(y)Q(y')) p = P (Q(y — y')Q(0)) p . 

In this article we have been mainly considering states which are in momentum eigen- 
states in the x coordinates, related to the y coordinates via (|2.2|). This does not lead to 
momentum eigenstates in the y coordinates. However, they do have definite momentum in 
the p- direction. To the extent that we can neglect other components of the momentum 
in the y coordinates we see that the charge correlator has a simple relation to the deep 
inelastic scattering amplitude and the ordinary parton distribution functions. In the gen- 
eral case we will need to evaluate expectation values of the form (p'\JJ\p). These require 



generalized parton distribution functions ||59|| . Thus, if we have a state with definite mo- 
mentum in the original x-coordinates, we will have a supersposition of momenta in the y 
coordinates and the charge two point function will be related to integrals over generalized 
parton distribution functions. We will not write a detailed expression here. 

Notice that the integral over x is divergent at small x. We think that this is due 
to the fact that the integral over y is also divergent for the charge correlator since the 
small angle singularity is not integrable. This divergence, though, is local in q and can 
probably be extracted without changing the overall picture. We have not checked this in 
detail. This problem is not present if we consider the energy correlation functions and the 
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relation to the deep inelastic amplitudes probed by gravitons. In that case all quantities 
are manifestly finite. 

The fact that in our problem we do not have ordinary plane wave wavefunctions in y 
has an interesting consequence. It was shown in ||] that, in the gravity regime, the leading 
power of q, which governs the short distance behavior in y, is controlled by a double trace 
operator. We will show below that this contribution is highly suppressed for operators 
that have definite momentum in x-space. We expect that this double trace contribution 
will also be suppressed at weak coupling when we consider plane wave states in x-space. 

Of course, everything we said here can be repeated for the energy correlation function, 
except that we should consider a deep inelastic process where we scatter gravitons from 
the field theory excitations. 



2.5. Energy correlations and the C parameter 

Let us make here a side comment on the relation between the energy correlators and 
other usually considered event shape variables. Event shape variables are certain functions 
of the four momenta of the observed particles which are infrared safe. One concrete example 
is the C parameter, defined as [^0| 



C=^ J d 2 n 1 d 2 n 2 £(n*i)£(n*2)sm 2 6 12 (2.50) 

where E is the total energy (and we assume that the total momentum vanishes). We see 
that the expectation value of C is given by an integral over the energy two point correlation 
function. 

On the other hand, it is common to compute the cross section as a function of C (see 
for example [pl[ ). This is just the probability of measuring various values of C, %. This 
calculation involves more input than the two point correlation function, since we would 
need to know all the moments of C, (C n ), to reconstruct 4^. 

The point of this short remark is to stress that, even though the C parameter is 
given by a product of energies, the computation of the cross-section as a function of C 
involves knowledge of the n point energy-correlation functions. Of course, in practice, 
is computed directly rather than going through the energy correlation functions. 
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3. Energy correlation functions in theories with gravity duals 

In this section we consider energy correlation functions in conformal field theories that 
have gravity duals. We first start with some general remarks on the energy correlators 
and the basic ingredients necessary to calculate them. Then, we will present explicit 
calculations for the energy one point functions, which are given in terms of three point 
functions in the gravitational theories. Finally, we add a general prescription for computing 
arbitrary n point functions. 



3.1. General remarks and basic ingredients of the calculation 

The general prescription for computing correlation functions of local operators in the 
CFT using the gravity dual was derived in |16|,|15| . Computations of the expectation values 



of the stress tensor for falling objects include @,^2[ . Since energy flux correlation functions 
are given in terms of stress tensor correlators, we simply need to perform the integral over 
time and take the limit in ( |1 . 1| ) . In order to simplify the computations it is useful to 
consider some coordinate changes. 

Let us start by writing AdS$ using the coordinates 

-(W' 1 ) 2 - (W ) 2 + (W 1 ) 2 + (W 2 ) 2 + (W 3 ) 2 + (W 4 ) 2 = -1 (3.1) 

The boundary of AdS$ corresponds to the region where W M — > oo. In that regime we 
can forget about the —1 in (|3.1| ) and we recover the coordinates Z M that we described 
around ( |2.3| ). It will be convenient to introduce three possible sets of coordinates which 
are natural from different points of view. The first two are 



Original: - = W' 1 + W 4 , = — , ^ = 0,1,2,3 

z z 

1 _ tj/O , w3 w-1 _ „.0 IJ/4 _ „,3 TXT". _ _ ^-2 



Easy : — = W° + W 3 , W' 1 = -y° , W 4 = -y 3 , W x 



(3.2) 



2 



2/5 2/5 

Of course the metrics are simply 

„ . . , , 2 dx 2 + dz 2 
Original : ds = ^ 

z 2 

Eag , ds 2 = -dy + dy~ + dyl + dy\ + dy 2 
asy. s ^ 
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(3.3) 



It is also convenient to introduce a third set of coordinates, which is defined as follows. 
We first choose three coordinates describing the H 3 subspace — (W ) 2 + (W 1 ) 2 + (W 2 ) 2 + 
(W 3 ) 2 = —r 2 for a fixed r 2 . The two other coordinates are chosen as 

W ± = W' 1 ± W 4 (3.4) 

Then r 2 = 1 — W + W~ and the metric is 

Hyperbolic: ds 2 = -dWUW~ - \ ^'^^^ + (1 - W+W~)d^ a (3.5) 

The advantage of this coordinate system is that it makes the SO(l, 3) symmetry of the 
problem manifest. This S'0(l,3) symmetry are the isometries of H3. In addition, the 
dilatation symmetry in the original coordinates becomes a boost in the plane, which 
is also a clear symmetry of the metric in this parametrization. 

The surface that is at the boundary of four dimensional Minkowski space can be 
extended to the interior in a unique way so that it is invariant under the symmetries that 
preserve the boundary of Minkowski space. In fact, this surface is simply given by W + = 0. 

The insertion of the stress tensor operator corresponds to a non-normalizable pertur- 
bation of the metric in the bulk. It will be convenient to derive first the expressions for 
the momentum in the y coordinates introduced in (|2.2|) . Since all the generators in ( |2.6| ) 
are given in terms of the integral of T--(y') over a line along y' , let us compute this first. 

I 12 

We insert T--(y ) on the boundary at y =0 and y = y =0 but at an arbitrary value 
of y' . We denote the boundary points with primes and the bulk points without primes. 
This induces the following fluctuation in the metric of AdS^, Qmn 9m n + ^mat, 

h MN dx M dx M ~ (dy+) 2 —- FrjSr- 2^ 2 T ■ u ( 3 ' 6 ) 

[-y+(y -y' ) + y( + y£ + y£ + ze] 4 

We can now perform the integral over y' . We use the formula 

1 1 

dy'~ - ,-^^. 14 ~<% + )73 ( 3 - 7 ) 
3 [y+y' +A + ieY A 6 

1 2 

Now, by a simple translation, we can set the energy operator at any other value of y' ' . 
We obtain 

h M NdX M dX N ~5(y + )(dy + ) 2 2 f 5 j^-xdy^dy^ = 

{yi + (yi - ViY + (2/2 - y 2 ) r 
3 

^5(W + )(dW + ) 2 T ^— j^— r^dy'M (3.8) 

{yi + \vi -Vi) + Km -V2) ) 

~S(W+)(dW+) 2 { ^ + ^ ) l ) dZ l dZ 2 
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where in the last line we have represented the boundary coordinates using (j2.4|) , with 
Z + = 0, in order to get the answer in a form that will make it easy to make coordinate 
changes. 

For example, if we wish to express the result in the hyperbolic coordinates ( |3.5| ) 
all we need to do is to express W and Z in terms of such coordinates. Let us define 
n = (ni, 77-2,71.3) to be a unit vector on a two sphere. On the surface W + = we can 
parametrize W° = cosh£, Wi = sinn^n^. It is then natural to take boundary coordinates 
Z° = 1, Z % = n' % . We take the limit ( — > 00 and we then have that 

W + W* ~ e<(l + n 3 ) ~ e< (Z" + Z>) , y[, 2 = = (3.9) 

The last equation is simply the change of coordinates ( |2.10| ). We then find the following 
expression for the generators 



E — > hf IN dX N dX M ~5(W + )(dW + ) 2 -— — — — -^dQ' 2 



1 



pi _^ h p MN dX N dX M ~ 6{W + )(dW + fn' l {wQ -±— dSY 2 (3.10) 

m -> h £ ^dx N dx M ~ 5(w + )(dw + r {w0 _^ w , )3 

In summary, we have computed the metric fluctuation that corresponds to the integrated 
insertion of the stress tensor that measures the energy deposited in the idealized calorime- 
ters that we are placing at infinity at some position n' on the two-sphere. In the original 
AdS coordinates the corresponding insertion would be localized on the horizon of AdS in 
Poincare coordinates {z = oo). We found it convenient to express the results in a couple 
of different coordinate systems that are regular at z = oo in order avoid having to take 
a limit. In these other coordinates we see that we are performing a measurement on the 
W + = surface. This amounts to sampling the wavefunction of the particles in the bulk 
at W + = 0. At W + = we have an H 3 subspace plus the null direction parametrized by 
W~ . The boundary of H 3 corresponds to the two sphere at infinity where we place the 
calorimeters. 

The form of the equations (|3.10|) does not make explicit the Lorentz covariance of the 
expressions. In order to see this explicitly we can rewrite them as 

P» - h% N dX N dX M ~ S(W+)(dW+r ( _^ z)3 d y (3.11) 
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where we sit at Z + = and we are integrating over a spacelike surface inside Y^=o Zfj-Z^ = 
0. The integration surface differential is defined to be such that dS u x Z x is parallel to Z v ' . 
Therefore the transformation of dS^Z x cancels the transformation of Z° . 

In the same way that we have discussed the graviton associated to energy flux mea- 
surements we can also consider the U(l) gauge field configurations associated to charge 
flow measurements on the boundary theory. The operator that corresponds to putting a 
counter at infinity at some specific location and measuring the charge corresponds to the 
following bulk gauge field configuration 

Q(n') ^A M dx M ~ dW+5(W+) w ^ WW2 (3.12) 

Having discussed the properties of the probe gravitons or gauge fields that represent 
our measurement, let us now turn to the field in the bulk that describes the state that 
we insert with the operator O. We can think of a scalar source to make things simpler, 
although our results will be quite general. We are interested in obtaining the field config- 
uration, <fi, in the bulk of AdS§ created by the insertion of the operator O of dimension 
A. If we insert the operator J d 4 x'4>o(x')0(x') on the boundary theory, then the bulk field 
configuration is given by [ |I6| , |1"5|1 



<p(x, z) = d x'4>o{x' 



z 



A 



[(x - x') 2 + Z 2 } A 
lim / d x'dnix') 



z> "o J " ■* ' (z') A (W.W') A (3.13) 



(W + = 0,W-,W) = / d 4 x 



D iq.x 



w_ _ W o x ,o + W i x ,i _ ie 



where we have first rewritten the result in a way that allow us to easily change coordinates. 
In the last line we wrote the expression for the bulk field at W + = in the case that 
4>q{x') = e iq ' x . Notice that we only need the wavefunction at W + = since that is where 
the graviton perturbation is localized. It is not hard to do the integral in (|3.13|) explicitly. 
There is, however, a very simple way to see what the answer should be. We are creating 
a state that is a momentum eigenstate. For the moment let us set q^ 1 = (q°,0). The 



13 This is completely analog to the fact in classical electrodynamics that the power radiated by 
an accelerating charge is a Lorentz scalar. 
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momentum generator corresponds to a bulk isometry generated by a Killing vector that 
becomes simpler at W + = 0, 

P£\w+=o = - 2iW »d w - (3.14) 

Of course, this is similar to the corresponding boundary statement ( |2.5| ). A wave function 
that diagonalizes all four of these operators has to be a plane wave in W~ and should be 
localized in the coordinates. In other words we have 

(j> q (W + = 0, W~, W) ~ {q°) A ~ 4 e iq ° w ~ /2 5 3 (W) (3.15) 

where W refers to a parametrization of the hyperboloid given by W l with i = 1, 2, 3. In 
these coordinates, W° is just a function of W since W^W^ = —1 at W + = 0. (q°) A 4 is 
a normalization constant that can be obtained from (|3.13| ) by considering the dilatation 
operator acting on both sides. In other words, 4>o(x') is dimensionless so that <fi scales as 
((Z°) A 4 - The overall constant in ( |3.15|) cancels out when we compute energy correlations. 
Note that in the y-AdS coordinates ( |3.2j ) the wavefunction with definite x-momentum, 
(|3.15|) , is localized at y\ = yi = 0, ys = 1 when y + = 0. 

Therefore, the general result is that an incoming plane wave (with no spatial momen- 
tum) gives us a very peculiar wavefunction that is S function localized in the H3 subspace 
at the origin W l = 0. In addition we find that the momentum in the W~ direction is 
proportional to the original energy. The wavefunction for an external operator with a 
generic value of the momentum can be obtained by performing a boost of this solution. 
The end result is again a wavefunction that is localized at a point i n H3. It is localized 
at = The momentum in the W~ direction is now _ V ^ 9 -. The wavefunctions 
corresponding to plane waves have a divergent norm since a plane wave wavefunction has a 
divergent norm. One can consider the regularized external wavefunction in (|2.28| ) . In that 
case we find a finite norm. We discuss this case in more detail in appendix B. As expected, 
one finds that the delta function is smeared over a region \W\ ~ — . We will continue to 
discuss wavefunctions for plane waves, but having in mind that we will eventually smear 
the 6 function in ( g3§ ), as in ( gggfl . 

Once we have the bulk wavefunction we can compute the energy flux one point func- 
tion. By considering the effects of the metric perturbation (|3.10[) , and considering an 
operator that creates the bulk wavefunction </>, we obtain the expression 



(S(n')) =N~ 2 [ dW-dZ 3 \ [(2id w -cf>*)(-2id w -mw^ 

J 4tt(W° - W.n'Y 

N 2 = [ dW~ [ dZ 3 {<P*(-id w -<P)+c.c. 



\w+=o 
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where S3 denotes the integral over the three dimensional hyperbolic space parametrized by 
W^, with WnW^ = — 1. The last factor, iV 2 , is simply the total production cross section 
and it is related to the two point function of the operator insertion. In other words, the 
two point function (OlC^OqlO) = ||O ? |0)|| 2 is the norm of the state. This norm is given in 
the bulk by the expression for N 2 . When we insert the wavefunction (|3.15| ) we see that a 
single point in the integral over hyperbolic space contributes. We finally get the expected 
result = O- 

3.2. Energy flux one point functions in theories with gravity duals 

Using the results above we are ready to calculate the energy one point functions 
for different type of sources. For a scalar source the symmetries imply the result ( |2.29|) . 
Because there are no free parameters we know this is the correct result and we did not 
need to go throught the previous discussion. The situation is more interesting for current 
sources. In a theory that has a gravity dual the three point function of two currents and 
a stress tensor can be computed from the bulk interaction between two bulk photons and 
a bulk graviton that follows from the bulk Maxwell action 

S = -J^2 J d 5 x^F 2 (3.17) 

where g is the bulk gauge coupling. This term in the action also determines the two point 
function. Thus, we can see that the three point function will be determined and we will 
get a particular value for in ( f2.30| ). We can find this value by noticing that for M = 4 
Super Yang Mills we had 02 = 0. Thus, any theory that has a gravity dual gives us a2 = 0, 
as long as the two derivative approximation ( |3.17| ) is valid. In general there will be higher 
derivative corrections to this action. Up to field redefinitions there is a unique higher order 
operator that can contribute to the three point function 

S = ~ J ofx^-gF 2 + J d 5 x^W^F^F Sp (3.18) 

where W tJ,u5p is the Weyl tensor. Here M* is some mass scale in the gravity theory 
determining the strength of the correction relative to the strength of the Maxwell term in 
the action. In order to see that this is the only operator that contributes we consider the 
possible three point vertices between two photons and a graviton in flat space. It turns out 
that there are only two possible structures. This onshell vertex is so constrained because 
there is no kinematic invariant that we can make purely with the external momenta, 
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which all square to zero. In fact the two possible interaction vertices consistent with gauge 
invariance are 

vi =e M „ [e?A£(e 2 .fci) + (1 <- 2) - K^ 2 (e 1 .e 2 )} 

(3.19) 

V2 =tp V k i k 2 ( £ i- fc 2) (ea-fei) 

where 2 are the polarization vectors of the gauge bosons and e Ml/ the polarization vector 
of the graviton. They are all transverse e\ ■ k\ = and = 0. The first arises from the 
quadratic action ( |3.17| ) and the second from the higher order correction in ( |3.18|) . 

We expect that the higher derivative corrections give us a deviation from a perfectly 
spherical energy distribution for the state created by the currents. Notice that the higher 
order correction will not contribute to the angle independent term in the energy one point 
function. The reason is that this term is related by the Ward identities to the two point 
function and the two point function is not corrected by the presence of the higher order 
operator in ( |3.18|) because the Weyl tensor vanishes is AdS. A detailed computation 
in appendix D shows that the higher order term indeed contributes to the anisotropic 
contribution to the energy correlation function 

Notice, in particular, that in non-supersymmetric weakly coupled QCD we expect 
that the higher derivative corrections are comparable to the radius of AdS since a 2 is of 
order one for weak coupling ( |2.32| ). 

This anisotropy is intimately related to the anisotropy in the gravitational field that 
is produced by a fast moving photon. Let us consider a photon with high momentum 
\p-\ ^> 1. We focus on the problem in flat space for the moment. Such a fast moving 
particle produces a metric of the form 

ds 2 = -dx + dx~ + dx 2 + 5(x-)(dx-) 2 h(x) , V 2 /i = (3.21) 

where V 2 is the flat laplacian in the transverse directions. This metric is an exact solution 
of Einstein's equations (with zero cosmological constant) and arbitrary higher derivative 
corrections [B31 . The particular form of the solution for h depends on the coupling of the 



photons to the graviton. For the lowest order action ( |3.17| ) we find that ho ~ ^ which is 
independent of the spin of the photon. Here we are focusing on the five dimensional case, 
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so that we have three transverse directions x. On the other hand, the second interaction 
( |3. 18|) gives a function h of the form 

». ~ = fe^t 11 ^ (3-22) 



\x\ 



where rij = Xj/|a:|. Note that this contribution to the gravitational field is sensitive to the 
spin of the photon and it has a quadrupole form. In the case that we have a large number of 
photons, this quadrupole tensor would be proportional to the polarization density matrix 
of the photons. 

Even though we've discussed the case of a photon, all that we have said so far can 
be extended to the case that we have a non-abelian gauge theory in the bulk, which 
corresponds to a non-abelian global symmetry in the boundary theory. 

We have a similar story in the case that the inserted external operator is the stress 
tensor itself. Then there are three possible vertices and three parameters specifying the 
stress tensor three point function. One of these parameters is fixed by the Ward identities 
and it multiplies the three point function that we expect from the gravity action. The 
other two parameters multiply higher order gravity corrections. In fact, the three possible 
gravity vertices in five dimensions are 

vi =^e\^ s e 3S p k p 2 + -/ ^ e\k{k p 2 + cyclic 

V2 =(k^\^ s k s 3 ) (e 3 p Xik%) + cyclic (3-23) 

Such vertices arise from terms in the action of the form 



J pi pi 



(3.24) 



This is one way to parametrize the higher derivative corrections. In principle we can have 
another curvature cubed term but it does not contribute to the three point function [|64 |. 



In fact, in |64[] such corrections were computed for various string theories. They found 
that 71 and 72 are non-zero in the bosonic string, only 71 is nonzero in the heterotic string 
and both 71 = 72 = in the type II superstrings. Incidentally, corrections to this action, 
yielding an R term, were computed for type IIB superstrings in AdS$ x S 5 in |55| and 



their effect on the 3 point function of stress energy tensors in M = 4 SYM was discussed. 
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One can compute the contributions of the higher order terms in the action to ti and 



t\ as defined in ( |2.37|) . After performing some calculation described in appendix D, we 
find 

t 2 48 ^i + / 72 v 

R AdS M pl R AdS M pl , v 

4320 72 



to leading order in the 7^. In addition we have assumed that the contribution of the 
operator to £2 will be smaller than the one from the W 2 operator. This is expected in 
the large radius limit because W 3 has more derivatives. For an Af = 1 supersymmetric 
theory £4 = 0. Using ( gjgg ) and (gjgg ) we get the expression for the R 2 coefficient that 



was derived in [EB. 



Notice that the presence of the first correction to the action (|3.24| ) can also change 
the angular independent part of the energy flux one point function. This change should 
be compensated precisely by a change in the stress tensor two point function in order to 
obey the Ward identity. 

We could also consider charge one point functions. In that case there are two (parity 
preserving) structures for the three point function [p|,^5|j67| . The coefficient of one of 
them is determined by the Ward identities and arises only when we have a non-abelian 
gauge symmetry in the bulk. It comes from the usual bulk term of the form J Tr[F 2 ]. The 
second structure arises from a bulk term of the form J Tr[F^ u F uS F s or more generally, 
from a bulk coupling of the form J f a b c F^ u F hV& F c s M with a totaly antisymmetric f a f, c . 
Notice that these terms do not necessarily come from a non-abelian gauge symmetry. They 
could come from a coupling between three different £7(1) gauge field strengths in the bulk. 

The parity odd terms come from Chern Simons couplings in five dimensions [|16| [j6T 



For example, for a gauge field we can have J A A F A F or its non-abelian generalization. 



3. 3. Comments on the n point functions 

After this discussion on one point functions, let us move on to n point functions. 
All we need to do is to consider metric fluctuations which contain several insertions of 
the energy flux operator. In general we would have to worry about the bulk tree level 
interactions among the bulk gravitons corresponding to the insertions of the operators. In 
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our case, there is an important simplification. This is due to the fact that the following 
plane wave solutions^ are exact solutions of Einstein's bulk equations [|69| , |70[ 



ds 2 = ds 2 AdS5 + (dW+) 2 5(W+)h(w) (3.26) 

where h(w) is a function defined on the transverse space, which in this case is a hyperbolic 
space H 3 of radius one, given by -(W ) 2 + (W 1 ) 2 + (W 2 ) 2 + (W 3 ) 2 = -1. The function 
h(w) obeys the Laplace equation on this hyperbolic space 

V 2 w h = 3h (3.27) 

Of course, one can check 

hH ' = 4n(W° - Wh^f (3 ' 28) 
is a solution and so will be an arbitrary superposition 

n 

h = ^h ri/] (3.29) 

3=1 

which represents the insertion of n calorimeters at angular positions given by n'-, j = 
l,---,n. This is summing all the gravity tree diagrams. We should now consider the 
propagation of the wavefunction on the background of this plane wave geometry. We want 
to consider the effects of each hf? to first order in its strength but the combined effect 

3 

of all of them. Let us recall how we would analyze this problem in fiat space first. We 
consider a flat space plane wave of the form 

ds 2 = -dx+dx~ + (dx+) 2 f(x + )h(x) + dx 2 (3.30) 

and we will eventually take the limit where f(x + ) — > S(x + ). The scalar field obeys the 
equation 

-4d-d+<f) - 4f(x + )h(x)d 2 _(p + V 2 </> - m 2 <p = (3.31) 

We now assume that f(x + ) is nonzero only within some small neighborhood of the origin 
— e < x + < e. This implies that f(x + ) varies rapidly. Thus we assume that this rate of 



14 For applications of this type of solutions in confining backgrounds see [68]. 

15 Here the normalization factor is fixed such that we obtain the total energy upon integration. 
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variation is much faster than the rate of variation of the wavefunction along the rest of the 
coordinates. In that region we can then approximately solve the wave equation (|3.31| ) as 



{x+ = e)=e-f-< f(x+)hd -<j ) (x+ = -e) 



[x 



-hd. 



(j)(x + = -e) (3.32) 



Generalizing this method to our case of interest we find that 



(f){W + = 6,W-,W») = e- hd ™-(t)(W + = -e,W-,W") 



(3.33) 



where denotes a point in H3. This is nothing else than a translation of magnitude h 
in the W~ direction. The same type of behavior was observed for scattering of particles 
off shock waves in JTT| and was used to study four point functions in the context of the 



AdS/CFT correspondence in [|72[| . 

The computation we want to do involves the overlap of the final state with the initial 
state in the background deformed by the insertion of the plane wave. In addition, we need 
to divide by the norm ( |3.16| ) . If we write h = J2j hfi'. an d we expand in each of the hfpj 
to first order we get the n point function 



(Heft)) =n- 



dW~dE:- 



w+=o 



(id w - </>* ) II H ( W ) i(- 9 w- T<t>] + cc. 

3=1 



(3.34) 



where the integral over dT,s is over the hyperboloid W^W^ = —1, \i = 0,1,2,3. Let us 
specialize this expression to the case that we have a plane wave external state, which leads 
to (gnSj) . In that case we find that all the hft, are evaluated at W — so that they 
become independent of the angle. Thus we get that not only the one point function is 
uniform but also all the n point functions are uniform as well. This implies that there are 
no fluctuations in the energy and an observer would see a uniform energy deposition in all 
the detectors. In other words, we get 



(em •••£«)> = (£)'' 



(3.35) 



This is what we would expect by thinking that fragmentation is very rapid at strong 
coupling as suggested in HID- For a state with a generic, but definite, momentum we find 



1 



(5 



2\2 



4tt (q° — q.n') 3 
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(3.36) 





which is simply the boosted version of the uniform distribution £ = that we get for the 
case where q = 0. 

This is the result for plane wave states. If one considers a generic state, then there can 
be fluctuations, but such fluctuations are parametrized by the fact that we have a wave- 
function for momentum. In other words, one can write the wavefunction (po(x) appearing 
in (|3.13|) in momentum space as <fio(p) = J d 4: xe~' ip ' x <J)o(x). We consider only wavefunc- 



tions which are nonvanishing in the forward light-cone p 2 < 0, p° > 0. We could consider 
other wavefunctions but the corresponding operators vanish when they act on the vacuum 
and will not contribute. Thus, in the formulas below we imagine that p^ is restricted to 
be in the forward light-cone. Then we can write the bulk wavefunction as 

(j)(W+ = 0,W-,W^)= -^X A e lXW ^MXW^) (3.37) 

Jo * 

We see that for a plane wave with purely timelike momentum we should set <fio = S(p° — 
q°)5 3 (p) and we recover ( |3.15| ). Inserting ( |3.37| ) into ( |3.34| ) we obtain 



IV n 10 A 

pip) =n- 2 ( p 2 ) a - 2 \Mp)\ 2 , n 2 = I A( P 2 ) A - 2 |0 i) i / ,)|^ 



(3.38) 



The factor of (p 2 ) A ~ 2 appears when we consider the norm of a state that has momentum 
p, see appendix A. This factor is determined by the dilatation operator. In appendix B we 
compute p for the localized wavefunction (po(x) given in ( |2.28|) . 



The final picture is that for a generic operator insertion we have a superposition of the 
results for each momentum, given by (|3.36|) with a probability weight given by p(p) which 



is giving us the probability of exciting the mode with momentum p^ 1 in the conformal field 
theory (or in the bulk gravity theory). 

For a generic (po(p) (|3.38|) gives non-trivial functions of the angles. The final picture 
for what we would see in each event is actually very simple. After we measure the energy on 
four of the calorimeters in each event, we can determine the value of p that is contributing 
and, therefore, the energies in all other calorimeters is determined. See appendix B for a 
longer discussion of this point. In other words, from event to event, we have some random 
variations which are completely captured by the distribution of momenta p(p)- In the bulk 
picture, we have a pointlike particle in the bulk with some wavefunction. Measuring all the 
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energies is tantamount to measuring the position of the particle on H% and its momentum 
in the direction W~ when it crosses W + = 0. We can view this as the horizon of AdS. We 
can say that we are simply measuring the momentum of the particle as it crosses the AdS 5 
horizon. In the approximation that we have a pointlike particle we have a small number 
of random (quantum) variables characterizing the event. We only have the position or 
momentum of the particle when it crosses the horizon. When we consider a string we 
have an infinite number of degrees of freedom and we can have much more variation in the 
energy deposition patterns. 

4. Stringy corrections 

In this section we study stringy corrections to the gravity results. First we consider a 
flat space problem that is closely related to the problem we need to solve in AdS. We then 
use these results to compute the leading order stringy corrections to the gravity results. 
Finally we study the small angle behavior of the two point function and we find the stringy 
version of the operator product expansion we discussed above. 




(a) (b) 

Fig. 4: (a) The AdS computation of the energy correlators involves gravitons 
that propagate from the boundary to the interior on an H3 subspace of the full 
AdSs space. The gravitons originate on the boundary of H3, at the point where 
the calorimeter is inserted, and propagate on H 3 to the interior, (b) Since the 
falling string state is localized on H3 we can approximate the computation by a 
flat space one. 

4-1. Strings probed by plane waves 

Let us first make the approximation that the AdS space is weakly curved and let us 
approximate the problem as that of strings in flat space, see fig. 4. In fact, we have seen 
that the state created by the operator insertion is localized on the transverse surface, the 
H 3 subspace. In addition, the energy flux operator corresponds to a graviton localized 
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at W + = 0. Thus, we can just look at the problem in a neighborhood of this point and 
approximate it as a flat space problem where we have a particle, or a string, with nonzero 
P- which it is being probed by gravitons with p- = that are localized along y + . Note 
that the probe gravitons are extended in the y~ direction. 

More explicitly, we can consider a flat space problem where we have a string with a 
non-zero value of p- which crosses a gravitational plane wave of the form 

ds 2 = -dy+dy- + (dy + ) 2 5(y+)h + df (4.1) 

where h is a function of the transverse coordinates obeying V 2 /i = 0. Due to the symmetries 
of the problem it is convenient to choose light cone gauge where y + = —2a'p-T. Recall 
that p- < is the momentum conjugate to y~ . Following the usual steps that lead to 
light cone quantization we find that we get the following light cone gauge Lagrangian for 
the transverse dimensions 

S= MJ dT J d ° d T[(dTy) 2 -(d a y) 2 ]-—p-J o dah(y(r = 0,a)) (4.2) 

Notice that the h dependent term is localized at r = 0, at a single value of the worldsheet 
time. Thus, the string propagates freely in flat space away from y + = 0, or r = 0. We will 
also assume that near y + ~ the string is localized near y = in the transverse directions. 
We can then compute correlation functions from the expression 

where |^) the full wavefunction of the string state in the light cone gauge theory at r = 0. 
We consider a function h which is a sum of a finite number of plane waves, h = hje lkjlJ ', 
and we expand to linear order in each perturbation hj. The mass shell condition is k 2 = 0. 
This implies that k is complex. (When we go back to the AdS problem it will be natural 
to take the component of k along the radial direction to be purely imaginary and the 
others to be real.) Let us assume that the string state corresponds to the ground state for 
the bosonic oscillators excitations, at least for the bosonic transverse directions where the 
momenta k are nonzero. We then find that we have to compute correlation functions of 
the form 



j j 

~ (-ip-r^cmiu^^ii I ^ni 2si 



(4-4) 
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where we have separated out the contribution from the center of mass and the oscilla- 
tors. Since the center of mass wavefunction is well localized we expect no contribution 
from it. Namely, we imagine a wavefunction which is localized near y = an thus we 
simply need to evaluate the plane waves in ( |4.4j ) at zero which just gives one. Namely, 
(0cm I Ylj e xkjy \ificm) ~ 1- Note that if we neglect the oscillator contributions we recover 
the gravity result following from ( |3.33|) . Therefore, the nontrivial contribution comes from 
the oscillators. Notice that these integrals are convergent if the /c's are all small enough. 
In the case of the two point function we have 

— |2sin-r hl - h2 = j= U * J = 
o 27T 1 2 1 v^F r(l + ^%^) (4 _ 5) 

In the second line, the 1 corresponds to the gravity result and the second term is the 
first correction. Naively, one might have expected the first correction to be of order a' . 
However, the first term vanishes and the order a' 2 term is the first non vanishing one. 

It is convenient to rewrite this result in position space. We find that the gravity result 
plus the leading order correction can be written as 



, r 2 

N~' z I dy~d 6 y 



y+=o 



i(j>*d°-(f> + c.c. 



hx(y)h 2 (y) + ^a' 2 (d i d j h 1 (y)d i dMy)) 



where N 2 = / dy d 3 y [i(p*d y -(j) + c.c] 

Jy+=0 



(4.6) 



where <fi is the wavefunction of the center of mass of the closed string state and h\ and h 2 
are two graviton plane wave states. We have also normalized the result. 

In fact, from our discussion we can easily see the origin of the a' corrections to the three 
graviton vertex in various string theories. These were computed in [64jj. We first consider 
the case where there is just one probe graviton in (^4.4|) . a' corrections can only arise if the 
initial state contains bosonic oscillators in the transverse directions. For a graviton in the 
superstring we have no bosonic oscillators in the initial state, only fermion zero modes. 
Thus the vertex is the same as the gravity one. In the case of the heterotic string the 
graviton contains fermion zero modes for the right movers and a bosonic oscillator for the 
left movers. Such an oscillator can give rise to momentum dependent terms of the from 
given by the second vertex in ( |3.23| ), but not like the third in ( |3.23| ). Finally, in the case 
of the bosonic string a graviton with indices in the transverse directions involves bosonic 
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oscillators for both left and right movers and gives rise to a vertex like the third in ( |3.23j ) 
(plus the first two, of course). 

It is interesting that the string result is finite. One might have worried that since we 
are using 5(y + ) wavefunctions we would obtain divergencies. As we will see in more detail 
below, the Regge behavior of the scattering amplitudes in string theory ensures that the 
results are finite. 

4-2. Leading order a' corrections to the two point function 

We now generalize this result to curved space. We can simply replace the ordinary 
derivatives in (|4.6| ) by covariant derivatives. However, in the AdS context h obeys an 
equation of the form V 2 /i = 3h, or more precisely V 2 /i = -^r — h, so that terms that 



AdS 



would have been zero in flat space are non-zero in AdS, so we seem to be faced with 
an ambiguity. However, these ambiguities only affect terms that do not have angular 
dependence, at least for the first correction. Thus, such terms only correct the constant 
part of the energy correlations. We can fix such corrections by demanding that we obey 
the energy conservation conditions. It is convenient to think about the problem in the 
hyperbolic coordinates ( |3.5| ). The graviton wavefunction associated to the insertion of a 
calorimeter at n' on the two sphere is given by (|3.16 ) 



1 1 

h ~ - ~ = (4.7) 

{W^-W.n'f (i + ™_#.^)3 

where we have have expanded the result around W ~ 0. We have already seen that the 
center of mass wavefunction is localized near the origin of hyperbolic space if we have a 
state created by an operator with zero spatial momentum on R > , see Q3.15 ). Thus, we 
can evaluate the derivatives in ( f4.6| ) and then set W = 0. This gives 



(S(n[)S(n' 2 )} = (|^) 2 



tv 2 a' 2 

1 + ~OA ~p4 \ d i d 3 hd i d 3 h \w=0 + COnst ) 
24 H AdS 



(4.1 



where the constant is an angle independent term that we cannot compute purely in flat 
space. It can be fixed so that we obey the energy conservation condition. In the end we 
find 

(4.9) 



{£{n' x )£{ti' 2 )) = (|V 



67T 2 . 2 „ 1 . 

1 + — (cos 2 12 -3) + 



for Af = 4 super Yang Mills, where cos 6*12 = n' x .n' 2 . We see that, as expected here, the 
distribution rises in the forward and backward regions. We have fixed the constant term in 
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the correction by demanding that the integral over one of the angles gives the total energy. 
The dots in (|4.8| ) denote higher order terms in the 1 / v^A expansion. 

In this derivation we have assumed that the state we are considering has no oscillators 
excited along the three transverse AdS directions. In the case of the superstring we can 
still have a massless mode with indices in the transverse AdS directions since those can 
be accounted for by fermion zero modes on the string worldsheet in light cone gauge. The 



result (|4.9| ) is very general and holds for any theory with a ten dimensional weakly coupled 
dual with an AdS§ factor if we replace 1/A — > a' 2 /R\ dS , under the assumption that we 
are creating a ten dimensional massless closed string with the external operator. 

4-3. Corrections to the n point function 

We now consider the n point function (£(ni) • • • E(n n )). We have seen that the gravity 
result is just a constant. Let us compute the stringy corrections. The leading deviation 
can be computed by expanding the full expression (|4.4|) up to quadratic order in products 
of ki.kj. The resulting correction is basically the same as the one contributing to the 
two point function ( f4.9| ). In order to see something new we can go to cubic order in the 
products ki.kj. In the end this gives us a correction to the n point function which looks 
like 



(£{n 1 )...£{n n )) = (j- 



(4.10) 







A 3 / 2 



[ ^ [nirij){njrt k ){ni.n k ) H ] + o(A 2 ) 



i<j<k 

where (3 is a numerical coefficientS and the dots denote terms that are necessary to ensure 
that the integral over each of the angles gives zero as well as a term that corrects the 
coefficient of the (n^nj) 2 term by an order A -3 / 2 amount. 

Thus, we find that for a strongly coupled field theory the energy distribution is uniform 
with small fluctuations which have an amplitude of order 1/vX In other words, 88 /S ~ 
-7=. The two point function of these fluctuations is given by the first non-constant term 
in ( f4.10| ). One might have thought that these flucutuations would be gaussian. However, 
we find that the three point function of the fluctuations is of order A -3 / 2 . Thus, when we 



16 



(3 = -1728 ^f- log[2 sin ^] l og [2 sin Si] log[2| sin <Si_^) y „ 51g ± 5 . 
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normalize the two point functions to one, the three point functions are of order one. Thus 
we conclude that the fluctuations are not even approximately gaussian. More explicitly, 
we can define a fluctuation operator 

5 = S ~Jf' , and 6 = y/XS (4.11) 

where the operator 8 is defined so that its two point function is independent of A. We then 
have 

(Hni)) =0 

(S(ni)S(n2)) =67v 2 [(n 1 .n 2 ) 2 + o(\~ 1/2 )} (4.12) 

(8(ni)6(n2)6(n 3 )) =P[(ni.n2)(n 1 .n 3 )(n2.n 3 ) H ] 

We see that the three point function is not parameterically suppressed relative to the two 
point function. Of course, they are both suppressed relative to the gravity result. 

4-4- Stringy corrections to charge two point functions 




(c) 

Fig. 5: (a) Feynman diagrams that lead to energy correlation functions. The 
gravitons do not interact before they touch the falling state, (p. We have also 
indicated the t and s channels as we define them in the text, (b) Diagram that 
leads to a divergence in the charge correlation function. The intermediate state is 
a graviton and the AAh vertex comes from the Maxwell term in the action. This 
divergence is cured by going to string theory and exploiting the Regge behavior of 
the amplitudes. In (c) we draw a diagram that can arise due to a higher derivative 
contact interaction in the gravity theory which could lead to a divergence in the 
gravity approximation. 

We now consider the two point function of a charge that is dual to a closed string 
mode. For example, we can pick one of the SO (6) currents in Af = 4 super Yang Mills. 
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More generally we consider a current associated to a symmetry that is carried by fields in 
the adjoint representation in the dual field theory. Imagine that the current comes from 
Kaluza Klein reduction. Then the corresponding vertex operator, in light cone gauge, has 
the form A+ — > d T ^e %k ' x where tp is one of the internal dimensions. We assume that the 
state we are measuring, does not have any oscillator excited in the <p direction and 
that it is not charged. Thus the one point function of the charge, is zero. The two point 
function is actually infinite in the gravity approximation. This is due to the Feynman 
diagram in fig. 5(b). This is intimately related to the fact that ( |3.12[ ) is not an exact 



solution of the gravity equations, but it sources a gravitational plane wave proportional 
to FL, which leads to the square of a S(y + ) function. We did not run into this problem 
in the gravity theory because there are no diagrams of this form due to the fact that the 
gravitational shock wave is an exact solution of the theory. Even in the gravity theory we 
could have run into this problem if we had had a higher derivative contact interaction that 
brings together two gravitons, as in fig. 5(c). 

Let us now compute the two point function in string theory. The corresponding flat 
space expression is similar to QL5|) , but with an extra factor coming from the contractions 
of the d T (f field coming from the two vertex operators. We get a result proportional to 



r 27T rlrr rr r >a 'k 1 .k 2 -2 TV— i -I- a ' fc l- fc 2 \ n > U U 

Jo (2tt)' 2 1 v^F r( o^) 4 



We have defined the integral by analytic continuation in k\.k2 We see that we get a perfectly 
finite answer in string theory. The string theory answer even goes to zero as we take the 
small momentum limit. Translating this flat space result fl4.13| ) to AdS as we did above 
we get 

(Q(ni)<2(n 2 )> = -^ni.n 2 = -^=cos# 12 (4.14) 

where 7 is a positive numerical coefficient. This result has the angular dependence that 
one would intuitively expect, with the two oppositely charged particles going in opposite 
directions. 

Let us emphasize once more an important point. Due to the fact that we are consid- 
ering shock waves which are highly localized, with 5(VF + ) wavefunctions, it is important 
to perform the computation in string theory rather than first taking the low energy limit 
of string theory and then doing the computation. We will revisit this point later. 
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A B 



A B 




x=0 



1>«- 



-XT=0 



(a) 



(b) 



(c) 



Fig. 6: (a) Worldsheet vertex operator insertions for charge correlators associated 
to closed string gauge fields. These charges are carried by the bulk of the world- 
sheet. In (b) and (c) we consider charges carried by open strings. We consider 
an open string stretching between two different branes called A and B. In (b) we 
consider the charge two point function for the U(1)a living on the brane A. The 
two vertex operators are inserted at the same point. In (c) we consider the charge 
two point function for charge living on the brane A, U(1)a, and the charge living 
on brane B, U(1)b- The vertex operators are inserted on different boundaries and 
the result is non-singular. 



Another interesting situation arises when we consider currents that act only on fields 
in the fundamental representation, such as flavor symmetries. In this case the currents 
live on D-branes in the bulk. For simplicity let us assume that we have two D-branes with 
two different U(l) gauge fields in the bulk. Let us call them U(1)a and U(1)b- We could 
imagine a QCD-like theory where U(1)a and U(l) b are different flavor number symmetries. 
At leading order in A we detect these charges in the detector only if we create a mesonic 
operator that contains the corresponding quarks. Consider a situation where we have a 
lorentz scalar meson where the quark is charged under U(1)a and the anti-quark under 
U(1)b- In such a situation we expect to find only one charge of each type in the detector. In 
this case the charge one point functions are (<2a(^)) = (Qb(^)) = j~- What is the charge 
two point function for U(1)a? From the boundary field theory point of view we expect it 
to be zero for generic angles since the charged quark can be detected only at one particular 
angle, since we are working to leading order in A where we do not create quark anti-quark 
pairs. On the other hand, to leading order in A, from the gravity plus maxwell theory 
in the bulk we get get a completely spherically symmetric distribution of charge. In this 
case, the divergent term coming from the Feynman diagram in fig. 5(b), is subleading in A 
and we do not consider it (string theory ought to make this 1/A correction finite too). In 
other words, the two point function for the charges is {Qa{^)Qa{^')) gravity = j^yz- This 
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contradicts the field theory expectations. The resolution is that the stringy corrections are 
so large that they completely change the gravity result. Let us first see how this works 
in the flat space case. Here, we can quantize the open string in light cone gauge and we 
will get an action very similar to the one we had for the closed string except that the 
photon vertex operator, which is inserted at r = 0, is also inserted at a = at one of the 
boundaries of the open string. We find that 

(^cm|e^ iy V^ 2?7 |^ cm )(0|e^ l2; - c(o ' o) e^ 22; - c(o ' o) |0) (4.15) 

If we ignore the oscillators we go back to the gravity result. However, the contribution 
from the oscillators involves a singularity, since both vertex operators are evaluated at the 
same point. Formally, this gives a contribution of the form 2a kl - k2 . If k±.k2 > 0, then 
we see that this vanishes. Thus, if we define the answer by analytic continuation we get 
zero for all values of k\.k<2,, including the physical values of k\.k2 for our problem (which 
are negative). The reason that stringy corrections have such a large effect is that we start 
with singular wavefunctions for the photon, which contain a 5(y + ). If we had started with 
a smooth wavefunction in the x + dimension we would have integrated the vertex operators 
along the r direction on the boundary of the open string worldsheet and we would have 
obtained a no n- vanishing function of k\.k^- 

On the other hand, if we compute the two point function for the two different U(l) 
charges, the charge carried by the quark and the charge carried by the antiquark, then we 
get the vertex operators at opposite points of the string and we obtain a finite answer 

(V'cm|e^ 1 ^ 2 ^|V'cm)(0|e^ 1 ^ c(0 ' 0) e^ 2 ^ c(0 ' fJ=7r) |0) ~ 2 2a ' klM (4.16) 



In this case the leading order a' correction to the two point function reads 



(Oi(ni)QB(n 2 )) = ' 



1 81o S 2 a 

1 -j=- COS 012 

v A 



(4.17) 



(4tt) 2 

We see that there is a tendency for the two charges to go in opposite directions, as one 
naively expects. Of course, at weak coupling the quark and the antiquark fly in opposite 
directions (if we have the simplest operator which contains only a quark anti-quark pair). 
If we were to consider a higher point function we would get zero again. 

The general lesson is that when we compute charge correlators it is very important to 
understand the effects of stringy corrections. 

Once we consider finite N corrections we do not expect the two point functions for 
the same U(l) to be exactly zero. 
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4-5. Small angle behavior of the two point function and the operator product expansion 

In this section we study the small angle behavior of the two point functions using 
string theory. The leading order correction to the energy flux two point function ( |4.9| ) is 
analytic at small angles, i.e. when #12 — > 0. As we will explain below this is no longer the 
case once all the a' corrections are included. We will show that at small angles there is a 
non-analytic term of the form |#i2| p with a power, p, that we will compute. This power is 
intimately related to the singularities in the first line of Q4.5Q as a function of k\.k2- 

Let us first understand how the singularities in the flat space answer ( [4.5| ) arise. These 
singularities are at o?k\.k% = —1 — 2n. We can rewrite this condition as 

2 + 4n 

t = -(fci +k 2 f = — , ra = 0,1,2,-.. (4.18) 

a' 

Similar looking singularities are a well known feature of string scattering amplitudes and 
they arise when an invariant, such as t, is equal to the mass of a string state. In that 
case we can view them as arising from the production of an on- shell closed string state. In 
our case, however, there are no states in the closed string spectrum with masses given by 
( rhl8|) . Thus we seem to have a puzzle. We will argue that we indeed have certain string 
states, but of a non-local kind. 




Fig. 7: (a) The poles of the ordinary closed string amplitude arise from the 
region where the two vertex operators are close to each other but are integrated 
over both r and a. (b) Wordsheet OPE for the problem we are considering where 
we have wavefunctions localized in x + . In light cone gauge this results in operators 
localized at r = 0. Thus we get singularities from the region of the integral where 
a 12 — * 0. 
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Let us first understand the worldsheet origin of these singularities. For concreteness, 
let us focus on the first singularity at o!k\.k 2 = — 1. We see that at this point the integral 
in (4.5) diverges at a = like J — . At a ~ the two closed string vertex operators of the 



external gravitons come close together, see fig. 7b. This looks similar to the ordinary OPE 
region of a closed string worldsheet which produces the usual closed string state poles. 
The crucial difference is that in our case the integral runs only over the sigma direction 
(see fig. 7b), while in the ordinary case it runs over r and a (see fig. 7a). For this reason 
the position of the poles has been shifted compared to the ordinary closed string poles. 
Schematically we have 

Usual Case: I dz 2 \z\ a ' klM 1 



a'h.k 2 + 2 (419) 

Our Case : / da\a\ a ' kl - k2 — 

J a'k\.k 2 + 1 

We have now understood how the singularities arise from the worldsheet computation. 
Before moving on, let us clarify further a confusing aspect of these singularities. All we are 
doing is to scatter four string states: the two gravitons, the state we are measuring and its 
complex conjugate. So, why are the singularities different than singularities of the ordinary 
closed string scattering amplitude? What happens is that we are choosing very peculiar 
wavefunctions for the two gravitons. These wavefunctions contain S(y + ) factors which 
implicitly carry an infinite amount of momentum. More precisely, in order to go from the 
usual momentum space result to our expression for S(y + ) wavefunctions we should integrate 
the momentum space result over fci+ and k 2 +. The four point amplitude is characterized 
by t = — {k\ + k 2 ) 2 and s = — (p + /ci) 2 , where p is the momentum of the incoming closed 
string state with non-zero p-. Since =0, we have that t is independent of ki+ and it 
continues to be given by the transverse components of ki, t = —(k\ + k 2 ) 2 . On the other 
hand s contains a contribution of the form s = 4p_/ci + + • • •. For the polarizations of the 
gravitons that we are choosing here the amplitude has the form |73|]0 

a 4 n , i\ r(i-^)r(i-^)r(i-^) 

A4 = P- - H — 7— — ; 77; — ; 7— , u = —s — t (4.20) 

V* Wr(i + ^)r(i + f)r(i + ^) ' v ; 



17 We take the momenta to be nonvanishing only in the first five dimensions. The falling string 
state is taken to be a graviton which has indices in the remaining five dimensions. Of course the 
two probe gravitons have indices in the ++ directions. 
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In the gravity limit we recover the results that we expect from the diagram in fig. 5(a) and 
a crossed version. We can take p- to be fixed. Then the integral over fcj + translates into 
an integral over s El. At large s, the four point amplitude is controlled by Regge behavior. 
The amplitude goes as 

A 4 ~ s~ 2+ ^ (4.21) 

So the integral over s converges at large values of s for small t 0. As we increase t, the 
integral over s first diverges when the amplitude goes like 1/s, A4 ~ 1/s. This condition 
is precisely the n = case in ( |4.18| ). We get the higher order ones by a similar reasoning 



by expanding the amplitude to higher orders in the 1/s expansion. One minor subtlety is 
that only even powers of 1/s give rise to 1/s terms in the full amplitude, after we adjust t 
appropriately!!. Odd powers of 1/s would lead to extra singularities beyond those given 
by ( |4.18[ ), see more on this below. Thus, we see that the poles ( [118] ) are associated to 
the high energy behavior of the string amplitude. This is to be expected since by probing 
the string at y + = we are taking a snapshot of the string state and this requires high 
energy scattering. The fact that the amplitudes we are computing are finite is related 
to the fact that the high energy scattering displays Regge behavior. In conclusion, the 
result we obtained in lightcone gauge is perfectly consistent with the usual structure of the 
Shapiro- Virasoro amplitude. 

A related remark that we can make at this point is the following. Let us go back 
to the case where we consider a neutral falling state probed by two closed string gauge 
bosons. Going to ten dimensions we can view the gauge bosons as Kaluza Klein gravitons. 
In that case the flat space amplitude is very similar to (|4.20[) except that we now have 



r(i - Y) r ( 1 - s r) r ( 1 "" 
r(i + ^)r(i + ^)r(i + ° " 



A = p- —, t-tI 77H K ( 4 - 22 ) 



If we take the small momentum limit of this amplitude we get a constant. If we integrate 
this constant with respect to s, in order to go to S function wavef unctions, then we get 



18 The integral over &2+ simply gets rid of the momentum conservation delta function in the 
fc+ direction. 

19 There are poles along the real s axis. As usual, we give these poles a small positive or 
negative imaginary part so that we are analyzing the amplitude in the physical sheet. Thus, the 
poles along the real s axis do not lead to divergences in the amplitude. 

20 In other words, A 4 ~ s - 2 +«'*/2 Cn {a't)/s n , but c 2k +i{a't = 4k + 4) = 0. 
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an infinity. This is the infinity that we mentioned above as coming from the fieid theory 
diagram in fig. 5(b). Of course the full string amplitude is not constant. Thus, once we 
go to string theory we should integrate the full string amplitude (|4.20|) , which goes as 
A g ~ s~ l+2 ^ L and converges if t is negative. We can define it by analytic continuation for 
other values of t. Thus, we see that for this particular case, taking the low energy limit of 
the amplitude first and then doing the s integral gives a very different answer than doing 
first the s integral and then the low energy limit, which gives (|4.13| ). In the case of energy 
correlations, if we first take the low energy limit of ( |4.21[ ) and then we do the integral we 
get the same answer as doing first the integral and then the low energy limit. In this case 
this happens because the contribution is coming mainly from the s ~ and u ~ region. 

It is useful to perform explicitly the worldsheet operator product expansion of the two 
graviton vertex operators, see figure fig. 7(b). We obtain0 

p_ e ife i-2/( T=0 ' cr )p_ e ife 2- 2/(0,0) ^ ^2 ^a'k 1 .k 2 ^ e i(ki+k 2 )y(0,0) _|_ 1 (4.23) 

The pole arises when the power of a is precisely 1/a. This gives rise to the n = case in 
( [4.18|) . The operator that appears at this point has the form 

p 2_ e ik. y m 2 = —k 2 = — (4.24) 

a 

This is the operator that appears on the worldsheet in light cone gauge. One is tempted 
to write an operator in conformal gauge that would reduce to (^4.24|) in light-cone gauge. 



Due to its peculiar p- dependence, we are forced to write an expression of the formi! 

(d a y + d a y + )U(y + y k -y (4.25) 

with k obeying the condition ( |4.24j ) . This operator is formally a Virasoro primary but 
is not a proper local operator on the string worldsheet. Similar operators were shown to 



control Regge physics in [74J] . Of course, our regime is closely connected with Regge physics 
so it is not a surprise that similar operators appear. The operator ( f4.25| ) , without the delta 
function, has spin j = 3 in the y + ,y~ plane, as we had in the field theory discussion. This 



21 In the following expression it is convenient to redefine the range of a to [— 7r,7r], such that 
insertions close to the operator at zero are given by small \o\. 



22 



Notice that the extra power of p_ appears to compensate the one appearing from the delta 



function 5(y + ) ~ 
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is related to the factor of p 2 _ that appears in ( 4.23 ). Notice that the complete operator 
( |4.25|) , including the delta function, has total spin 2. It corresponds to the field theory 
operator Uj-i with j = 3. 

The operator (|4.24|) is the leading contribution in the worldsheet OPE ( }4.23| ). As 
we expand the exponentials at higher orders we pick up new operators which contain 
derivatives with respect to the transverse directions. Some of these operators can have 
transverse spin. These strings states have higher masses. They all have spin j — 3 in the 
y directions. Notice, however, that terms that come with odd powers of a in the ... in 
( ET23T ) vanish upon integration over a. The reason is a Z2 symmetry from interchanging 
a — > —a. This is completely analogous to the fact that terms that are not symmetric under 
the interchange of z —>■ z in the usual case where we integrate over the whole complex plane 
vanish upon integration. This is nothing else than the level matching condition. These 
terms are the same as the ones discussed above in connection with the singularities for odd 
powers of i. Now we understand why these poles are absent from ( f4.5| ). 

In the end, the singularities arise from a fairly ordinary worldsheet operator product 
expansion in light cone gauge. On the other hand, we cannot associate the corresponding 
worldsheet operator to an ordinary closed string state. This is related to the fact that 
the operator product expansion of two energy flux operators in the field theory lead to 
non-local operators in the y~ direction. In the field theory we were not too disturbed by 
the appearance of operators that are non-local in the y~ direction, so we should also not 
be surprised that in string theory we also get string states that are non-local in the y~ 
direction. These string states are localized at y + = 0, carry zero p- and are local in the 
transverse directions. 

The final conclusion of this discussion is that we should interpret the singularities in 
( |I.5|) as arising from the propagation, in the transverse space, of non-local string states 
created by the operators like ( |4.24J) or ( ]4.25| ). 

We will now argue that the short distance singularity of the energy flux two point 
function is governed by the operator associated to the first singularity in (|4.5|) . Since we 
are working in the regime of large but fixed A we might imagine that we could always 
expand the two point function as in the second line of ( |4.5[ ) . This is not correct if the angle 
is very small. In that case the relevant relative momenta are of order t ~ la 1 12 ~bt — • Thus, 

|012| H AdS 

at angles of order 6 ~ A -1 / 4 we cannot use the approximations leading to fl4.9| ). However, 
we can use the interpretation given above to the poles in t to write the flat space result 
in the first line of ( |4.5|) as a sum over contributions of poles. Then each pole corresponds 
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Fig. 8: Small angle expansion of the energy correlation function. The expansion 
is dominated by the propagation of a spin three non-local string state, denoted 
here by a thick red line. 



to the contribution of a physical (but non-local) string state that is localized in the y + 
direction but propagating in the transverse directions. We generalize this result to AdS 
by replacing the transverse space by H 3 . Now the non-local string states propagate on 
the H 3 subspace of AdS 5. These states propagate from the center of H 3 , where the string 
state created by the localized operator insertion is concentrated, to a region near the H3 
boundary, near the insertion of the two energy flux operators, see fig. 8. At large distances 
from the H3 center we expect that the wavefunction of the non-local string state goes as 
1/\W\ A , \W\ > 1, with 

A ~ mR A dS ~ v^A 1 / 4 + ■ • ■ (4.26) 

where m is given, to a good approximation, by the mass of the flat space state computed 
in (|4T24D . Incidentally, we can calculate the conformal weight A of other (generally non 
local) operators with arbitrary spin in the same manner. They correspond to the string 
states 

(d a y + d a y + )i5(y+y k -y (4.27) 
The mass of these states is given in flat space by m 2 = —k 2 = \{j — 2). Therefore 

A(j') ~ V2 Vj^2A 1/4 + ■ ■ • (4.28) 

This formula is expected to be a good approximation only for j <C A 1//2 , since it was derived 
assuming the flat space approximation. For very large values of j we get a logarithmic 
behavior in j, see [75]. Of course this is simply the analytic continuation of the leading 



Regge trajectory. 
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The dots in Q4.26Q and ( |4.28|) denote terms independent of A as well as higher order 
corrections. We then need to compute the overlap of a wavefunction which decays like 
1/|W| A for large \ W\ with the wavefunctions of the two gravitons associated to the energy 
flux insertions. We find a behavior 

(£(6 1 )£(6 2 ) ■ ■ ■) ~ ^ 2 - 6 (W 3 - 1 (^ 2 ) • • •) (4.29) 

where U3-1 is related to the lightest spin 3 non-local operator, with zero p_, which at 
strong coupling has a large dimension ( }4.26|) . In string theory this expectation value is 
computed by inserting the operator ( |L24|) . 

In conclusion, the structure of the OPE is precisely what we expected from general 
principles in any conformal field theory. At weak coupling the operator product expansion 
is dominated by operators of twist slightly bigger than two. This leads to correlation 
functions that are highly localized along certain jet directions. For any value of the coupling 
the operator, or string state, that dominates has zero p- and spin j = 3. At strong 
coupling, the operator acquires a large twist given in ( f4.26|) . The fact that operators with 
spin j > 2 have large dimensions at strong coupling is seen to be intimately related with 
the fact that the energy distribution is uniform. Of course, this fact is also connected with 
the validity of the gravity approximation in the bulk. 

5. Summary, Conclussions and open problems 

Let us summarize some of our results. 

We studied energy correlation functions in conformal field theories. Energy correlation 
functions are an infrared finite quantity that is useful for characterizing the states produced 
by localized operator insertions in a field theory |)7|,fl0lJl3|l . 

They can be computed for all values of the coupling since they involve the stress 
tensor operator and make no reference to a partonic description. This is more manifest 
at strong coupling where the partons are difficult to see in the gravity or string description. 

After a conformal transformation these energy correlation functions amount to mea- 
suring the state along a null surface. More precisely, each "calorimeter" insertion corre- 
sponds to an integral of the stress tensor along a lightlike line, f dy~T--, ( |2.6|) . 

We have argued that the small angle behavior of the energy correlation functions is 
controlled by an operator product expansion which features non-local light-ray operators 
of definite spin. When two calorimeters come close to each other we have a spin three 



55 



operator {£ (9i)£ (#2) ■ ■ •) ~ |^i2 1^ 3— 4 {^3—1 (6*2) • • •)• These operators can be discussed for 
any coupling. We recalled the weak coupling expression for the twist [|33| (|2.22| ), and we 
also computed the twist at strong coupling r ~ v^A 1 / 4 O), after having identified the 
string states that are dual to the operator U3-1. These are not ordinary closed string 
states. They are peculiar string states localized along x + = that have non-local vertex 
operators on the covariant worldsheet but do have a local description on the worldsheet 



in light-cone gauge. Closely related string states appear in the Regge limit WM - Despite 
their unfamiliar features they control the short distance singularities of energy correlation 
functions. 

The light-ray operators that appear in the small angle behavior of the correlator are 
related to the ones that control the moments of the parton distribution functions. In 
fact, one can write a precise relation between the energy correlation functions on a special 
state and a particular moment of the functions that govern the deep inelastic scattering 
amplitude (|2~49l ). 

We have seen that energy flux one point functions in states created by currents or 
stress tensor insertions have an "antenna" patterni! which is determined by the three 
point functions in the conformal field theory ( |2.33| )( |2~3"T|) . In the gravity description this 
pattern is spherical but as we include higher order corrections to the gravity action we start 
seeing deviations from the spherical pattern ( |3.2U| ) ( |3T2l)| ). These deviations are sensitive to 
the spin of the operator that created the excitation in the conformal field theory. In the 
particular case of M = 4 super Yang Mills, the energy one point function is spherical for 
all values of the coupling. In more general Af = 1 superconformal theories we find that the 
antenna pattern, (|2.36|) ( |2~39| ), depends on the parameters a and c that characterize the 
three point functions in the current/stress tensor multiplet ||48|1 [ I7| . These results are exact 
expressions, valid for any coupling. They depend only on the two anomaly coefficients a 
and c defined in [58]. Demanding that the energy that calorimeters measure is positive 
we get a constraint on a and c, \a — c\ < c/2, which is saturated for free field theories (of 
course, c > 0). 

We gave a general prescription for computing the energy correlation functions on the 
gravity side. The operator insertion in the field theory produces a string state that falls 
into the AdS horizon. Energy correlation functions depend on the wavefunction of this 



23 Remember that individual events do not present this pattern. This refers to the one point 
functions which consist of averages over events. 
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string state at the AdS horizon. The falling string is probed by particular shock waves 
associated to the insertion of each calorimeter. This can be computed in a simple way by 
choosing a coordinate system in AdS$ that is non-singular at the horizon. We can view 
the computation of the energy correlation functions as taking a snapshot of the falling 
string state as it crosses the horizon. In the gravity approximation the result depends 
only on the momentum distribution of the initial state and it is independent of the spin 
or any other property of the string state we consider. If the state carries a purely timelike 
momentum = (q , 0), then the energy distribution on the detector is perfectly spherical 
with no fluctuations. As we include stringy corrections we find small fluctuations that 
are inversely proportional to the square of the radius of AdS in string units (or 1/va) 
([OD( |4.lOD . These fluctuations are small but they are not gaussian ( |4.12|) . Since the 
shock waves we are considering are infinitely localized one might worry that this leads to 
divergences. In fact, they would lead to divergent answers in a field theory context (at 
least in some cases). The Regge behavior of string amplitudes at large energies ensures 
that the results we obtain are finite. 

It should be fairly straightforward to generalize this discussion to other dimensions. 
The discussion in 2+1 dimensions might have some condensed matter applications, similar 



to irai 



It would also be interesting to understand finite N corrections. 

There has been a great deal of progress in computing perturbative scattering ampli- 



tudes in M = 4 super Yang Mills, see [77] for example. From these scattering amplitudes 
one can compute the energy correlation functions. On the other hand, since the energy 
correlation functions are already infrared finite, it would be nice to see if any of the meth- 
ods developed to compute amplitudes could be extended to compute the energy correlation 
functions directly, without having to compute the amplitudes first. 

Another possible direction would be to consider the "hadronization" corrections for 
a non-conformal theory. In particular one could imagine a non-conformal theory with a 
gravity dual. In a confining theory with a gravity dual we expect that these corrections 
will be large in the large N limit because the strings cannot break. One could also try to 
understand situations where the theory becomes free in the IR, such as 4 + 1 dimensional 
super Yang-Mills, or dimensionally regularized M = 4 super Yang Mills. 

It would also be interesting to generalize this discussion to more complicated initial 
states such as the one resulting from the collision of two closed string modes in the bulk. 
This would be analogous to pp collisions. 
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Finally, this discussion might have some implication for black holes, since energy 
correlations are a way of measuring the final state of Hawking radiation and its non- 
thermal properties. 
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Appendix A. Positivity of J dy T__ 

Let us consider first free field theories. The classical expression for the stress tensor for 
the Maxwell field, T__ ~ Y2i=i 2 * s explicitly positive since it is a sum of squares. 

On the other hand, the quantum expectation value of T__ can be negative. Let us recall 
why this happens. Formally, we also have the sum of squares of hermitian operators, so 
that we would also expect a positive answer. However, when we normal order we subtract 
and infinite constant. Then the normal ordered expression is not a sum of squares of 
hermitian operators. In fact, we have schematically T__ ~ (o^) 2 + a^a + a 2 where we 
have separated the operator into terms with different numbers of creation and annihilation 
operators. By considering a state of the rough form = |0) + ta[a\ |0), and using that 
the vacuum expectation value of T__ is zero we find that (^|T__|^) ~ i?e[ce] + o(e 2 ) 
where c is some number. By taking e to be a small complex number we see that we can 
make T_ _ negative at a point |78[ . 



Let us now consider the integrated expression £ = f dy This expression has 

the schematic form 

/■OO 

E~ / dp+p+(a p+ (y)ya p+ (y) (A.l) 
Jo 

we thus see that we have the integral of products of operators and their adjoints. This is an 
explicitly positive operator. We have used the variable p + ~ — 2p- which is positive. Notice 
that terms with two or two a operators have disappeared from ( |A.1|) due to the following 
argument. The integral over y~ enforces that the total p + should be zero. However, 
creation operators can only increase p + , thus we do not obey the p + = constraint with 
only creation operators. Further discussion on the null energy condition for free fields can 
be found in |L9| . 



Let us consider now an interacting field theory. If we choose the gauge A_ = 0, 
then the stress tensor T__ continues to be quadratic in the fields and the above argument 
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would hold. Of course, this argument is not too convincing since we might be ignoring 
renormalization subtleties or problems with the gauge choice. It would be nice to find a 
more general and solid argument for an interacting field theory. 



Appendix B. Energy distributions in gravity for generic states 

B.l. Energy distributions for general states 

Here we show how to go between a discussion of n point functions and the computation 
of probabilities for seing various energy distributions on the detector. Sometimes one might 
be interested in computing the probability functional p[£(9)] for measuring a particular 
pattern of energy deposition on the calorimeters. When one computes jet amplitudes one 
is computing probabilities of this kind, where one integrates over certain regions, such as 
the low energy region between two jets, etc. 

If we are given p we can compute the n point functions, as (£(9i) ■ ■ ■ £(8 n )) = 
JV£p[£]£(8 1 )---£(9 n ). 

Formally one can also go in the other direction by computing the generating functional 
for energy correlation functions, (e % I d e W) £ ( e )y This expression is a functional of X(6) 
and its expansion in powers of A gives us the n point functions. Then p is given by 

p[£'\ = J V Xe^I d20x ^ £ '^(e l I d2em£ie) ) (B.l) 

Just in order to see how this works, let us start with the n point functions given in 
( |3.38| ). We can easily compute the expression 

(e< / = J d 4 qpiq)e ifd>ex(e)E qf >(0) (R2) 



where £ q i*(9) is the function in (|3.36|) and p is defined in ( |3.38|) . After doing the functional 
integral in ( |Ij.l| ) we get 

p[£] = I d 4 qp(q)l[8[£(d) - £ q (9)] (B.3) 
J e 

We see that we have a continuum of S functions, one for each angle. But we are integrating 
only over four variables. Thus, once we fix the energy at four points, the energy at all 
other points is also fixed. 

The general, formal, string theory expression for the energy correlators has a similar 
form, except that we have to integrate over the infinite number of variables specifying the 
string wavefunction. At finite N we would also have many string states. 
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o (t 2 +S 2 ) 



„0. {s a „-„0.t (t 2 +S 2 ) 



B.2. Bulk wavefunction for a localized state 

Here we start with the wavefuntion 4>o(x) ~ e~ iq t e ^ that we mentioned in 
( |2.28|) . Its Fourier transform is 

O (p) = / d 4 xe- ipx M%) = I d 4 xe ipVt - l P s e- iqUt - 

J J (B.4) 

^4 e -^[(/- g V + ( P f] 

The bulk wavefunction then has the form 

rOO , q-2 

<P(W+ = 0, W~,Wn ~ (q°) a / rfAA^ 1 ^^ /^-^^ - 1 )^!^! 2 ] (B.5) 

./o 

We are considering the case that q°a ^> 1. We then see that as soon as \W\ ^> l/(aq°) 
the answer is exponentially suppressed. In the region with large \W\ ~ W° > 1 we can 
do the integral QB.5Q by saddle point approximation and we find 

<j>(W) ~ (q°) A (W°)- A j+jrr e^-/(2H") e -i# (B . 6) 

(aq ) 

We then insert this in (|3.34| ) to find an expression of the approximate form 

(£(n;)£(n' 2 )) ~ gge" * / dS 3 /T „ mA , 9 * ^ (B.7) 

where we used that iV 2 in ( |3.34|) is not exponentially small, and in (p.7|) we have kept only 
the leading exponential behavior in aq°. We have also approximated the integrand in the 
large W° region which we expect to dominate for the singular small angle behavior of the 
two point function. Finally we find the singular small angle behavior 

(n 2 )) ~ |#|2A-4 e - ((7g ) 2 /4 (Rg) 

This is precisely the power we expect for the double trace contribution as was discussed 
79fl . This term is exponentially suppressed when we consider a state with definite mo- 



rn 



mentum. Thus, the term that gave the largest contribution in the deep inelastic scattering 
analysis in W§1 does not contribute to energy correlators when we consider states created 
with definite momentum. They do contribute if the state does not have definite momen- 
tum in x-space. In fact, we saw that a state with definite momentum in y space is directly 
connected with the deep inelastic scattering amplitudes ( |2.49| ) . Such a state does not have 
definite momentum in a;-space and will not have the exponential suppression that we get 

in god. 



We should emphasize that the contribution to (|B.8|) is coming from the region where 
the particle is crossing the horizon at a position that is close to where the calorimeters are 
inserted. 
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Appendix C. Computation of the energy and charge one point function in 
terms of the three point functions of the CFT 

We denote by O any operator, which could be a scalar operator S or a vector e.j or 
a tensor tijTij. 

Let us start by recalling some formulas for two point functions 

1 



(o\s(t, x)S(o, o) |o> =— — — — -a 

[—{t — l€) z + \x\ z \ 

(0\T(S(t, x)S(0, 0))|0) = - \ , iA 

[-t 2 + \x\ 2 +ie] 



(CI) 



where the first is not time ordered and the second is time ordered. Of course the same 
prescription works for vector or tensor operators. The operator insertion with a definite 
timelike momentum q^ = (q°, 0) can be written as 



ego) = J dt e - iq0t o(t)\o) (c.2) 



and it creates a state with energy E = q° > 0. The fourier transform of the two point 
function is 



/ 



d 4 xe~^ x x = c(A)9(q°)(-q 2 ) A - 2 , c(A) = * } (C.3) 



This is the norm of the state that (|C.2| ) creates. This will also give us the total production 
cross section if the operator O couples to the standard model. As remarked in |8(J the 
positive norm condition implies A > 1. 

We are interested in starting from the ordinary expressions for the correlation func- 
tions in position space and extracting the limit that corresponds to the energy or charge 
correlators. In doing so, it is important to order the operators appropriately. For the 
non-time-ordered three point function the correct prescription is 

(0\S(x 2 )S(x 1 )S(x 3 )\0) = 

= 1 (C4) 

{[-(fas - + (x 23 ) 2 ][-(ti3 - ^ e ) 2 + (xi 3 ) 2 ][-(t 2 i - *e) 2 + (x 21 ) 2 ]} A / 2 

If one considers tensor operators we get similar denominators and we choose the same 
ie prescription. This ie prescription is a simple way to enforce the right ordering of the 
operators. Another way to say this is that an operator that is to the 'left' of another should 
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have a more negative imaginary part in the time direction. When one does perturbation 
theory, it might be convenient to use time ordering along a Keldysh contour. However, for 
our purposes this simple prescription suffices. 

Let us first show how to extract the energy correlation for a state created from a 
scalar operator with fixed momentum, or at least fairly well defined momentum, as in 
( j2.28| ). In this case, we know that the answer is independent of the angles and that the 
overall coefficient is determined by energy conservation. Nevertheless it is instructive to 
discuss this case in detail since the computation is the simplest and one can apply a similar 
method for other cases. Our method is not too elegant and there is probably a more direct 
and elegant method than the one we applied here. 

We extract the energy correlation by directly performing the limit and the integral 
in ( |1 . 1| ) . We use translation invariance to fix the position of the first operator at xs = 0. 
For simplicity we place the detector along the direction z, so that x\ = (£, 0, 0,r). We 
will take the limit r — > oo. If t is generic, then the three point function will decay as 
1/r 8 since it would be determined by the dimension of the stress tensor and the operator 
product expansion. This would be decaying too rapidly in order to give a finite large r 
limit. However, there is a larger contribution from the region t ~ r, the region on the 
light-cone of the inserted operators. This is the region that will contribute. Of course, 
this is precisely what we would expect in a theory of massless particles. It is convenient to 
define coordinates x = t ± r. We will find that the region with finite x~ will contribute. 
In addition, only the T__ component of the stress tensor can contribute. The integral over 
t can be traded for an integral over x~ . So we first take the r —>■ oo limit. We then do the 
integral over x~ and at the end we do the integral over x^. 

Let us see what each of these steps gives us. In order to follow this appendix the 
reader would need to have a copy of the paper by Osborn and Petkos HH], since we will 
make frequent reference to it. We start with the correlation function for 

(0|5(x 2 )T_.( a ; 1 )5(x3)|0> ~ 2A J 2 2 ( ^ - ^) (C.5) 

x 23 ^12^13 V X 12 X 13/ 

from equation (3.1) of |ff5 |. We are not going to keep track of overall numerical coefficients. 
After multiplying by r 2 and taking the r —>■ oo limit we get 

^<°l W-MSMIO) ~ ^ J — w J—— f (C.6) 



r 
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We now perform the integral over x . Note that we can close the contour on either the 
upper or lower x~ plane and pick one of the two poles in ( |C.6|) . We then find 

hm r 2 (0|5(x 2 ) J rfx"T__ ( Xl )S(0) |0> ~ ( C " 7 ) 



We now integrate over the two transverse x<i coordinates and use the wavefunction in ( C.2|) 
does not depend on them. We find 



lim r 2 / d 4 x 2 e iq ° t2 (0\S(x 2 ) [ dx"T__(xi)«S(0)|0) 

r^oo J J 



I 



dt 2 dz 2 e i ^ r - 1 2 - 3_— 3 ~ %°)(3°) 2A - 3 

[-(i 2 - + z|] A 2 (t 2 - -22 - 2ze) 3 

(C.8) 

where we have also done the remaining two integrals. When we divide by the two point 
function ( |C.3| ) we get {£) ~ q°. The numerical coefficient can also be computed at each 
step. Of course, this gives the right answer (£) = |- due to the Ward identity which fixes 
the coefficient of the three point function ( |(J.5| ) in terms of the coefficient of the two point 
function (|Cp , see eqns (6.15), (6.20) of 

This procedure can be repeated replacing the operator S by a current e.j. The com- 
putations are identical but with more indices and we use a computer. In this case the three 
point function of a stress tensor and two currents is fixed by conformal invariance, plus the 
Ward identity, up to one unknown coefficient. Conformal invariance leaves two possible 
structures and the Ward identity fixes the coefficient of one of them. In this case we find, 
as expected, that the energy correlation function depends on the angle with respect to the 
vector e. Here we simply quote the value of the parameter a 2 , introduced in ( ^.3U| ), in 



terms of the parameters e,c defined in (3.13) and (3.14) of [5S]c3. We find that 



_ 3(8e - c) ^ 3 EMLzM!L (CQ) 

2 2 (e + c) EM) 2 + ^f) 2 

where we indicated the value for a free theory with bosons and Weyl fermions of charges 
and q™^ . The combination (e + c) is fixed in terms of the coefficient of the two point 
function of two currents via 6.26 of 114511 . 
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Here we added hats to the parameters e and c in [45] so that they are not confused with 



other parameters in the present paper. Note also that [45] uses some of these letters with multiple 
meanings through their paper. 
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We can do this also for the correlation functions of the form {0\e*jTij£eijTij\0). We 
can then compute the coefficients t 2 , £4 introduced in ( |2.37| ) 



_30(13a + 4o-3c) 15(-4n v +n wf ) 

14a - 2b - 5c (n b + 12n v + 3n wf ) 

15(81a + 326 - 20c) 15(n b + 2n v - 2n wf ) 

2(14a-2o-5c) 2(n b + 12n v + 3n wf ) 



(CIO) 



where a, 6, c are defined in (3.19)-(3.21) of PB| . We have also indicated the result for n v , 
n b and free vectors, real bosons, and Weyl fermions (one complex dirac fermion would 
give n w f = 2)11. Again, the combination appearing in the denominator is fixed in terms 
of the stress tensor two point function, see (6.42) of [fl5 . 



Two combinations of these three coefficients are related to the values of a and c defined 
through the conformal anomalyil 

Tfi = -^W 2 - -^E (C.ll) 
16tt 2 16tt 2 V ! 

where W is the Weyl tensor and E = R^s P R fJluSp - AR^ V R^ V + R 2 is the Euler density, c 
also sets the two point function of the stress tensor. The coefficient a can be expressed in 
terms of the three parameters in ( |(J.1U ) as 



a _ (9a -2b- 10c) _ 2n b + 124n„ + lln wf . . 

c 3(14a-26-5c) 6n b + 72n v + 18n wf 



This follows from (8.37) in |(45[| . The values for a free theories were computed in [^5, jQj. 
From the positivity conditions ( |2.38|) it is possible to get general bounds on this ratio. We 
find 

H > a - > I (C.13) 
18 - c - 3 V ; 

where the lower bound is saturated by a free theory with only scalar bosons and the upper 

bound by a free theory with only vectors. Notice that this bound holds for any conformal 

theory while the more restrictive bound ( 2.40|) holds for supersymmetric theories. We can 



also add that for a M = 2 supersymmetric theory we can find a similar bound by using 



25 d, b, c here should not be confused with the e, c of the previous paragraph. In particular the 
two c are not the same 1143 1. 



20 



Do not confuse these a and c with the parameters a and c of the previous paragraph. 
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as limiting cases free theories with only vector supermultiplets (- = |) and free theories 
with hypermultiplets (~ = ^). Therefore | > ~ > |. This agrees with results from [|54| . 



In an A/" = 1 supersymmetric theory there is a relation between the three parameters 
d, 6, c which is obtained by setting t 4 = in ( |C.10|) . In this case, the two coefficients 



in (|C.11| ) specify completely the three point functions of the stress tensor. In a non- 



supersymmetric theory, we have one more parameter beyond the two in ( |C.ll 



Finally, we can repeat this exercise for the correlation function of three currents to 
find that the coefficient introduced in ( |2.43|) 



3 5a-4fr ^ ^C(n) b -C( ri ) wf 
2 2 (a + 46) £ i C(r i )6 + 2C(r^/ ^ ■ > 



where d, b are defined in eqn. (3.9) of ||45| , and are not the same as the ones in the 



previous paragraphs. Here are the representations of the bosons and Weyl fermions. 
And C(r) is defined as tr r [T a T b ] = C(r)5 ab . Again, the combination (d + 46) sets the two 
point function of the current. And a-i vanishes in a supersymmetric theory since there is 
only one structure contributing for the supersymmetric case [^J and it vanishes for a free 
supersymmetric theory. 

One can also take similar limits of the parity odd part of the three point function of 



three currents and one obtains the result in (|2.44j). Finally, starting from the correlation 



functions for two stress tensors and a current derived in |56|] one can derive the charge 
distribution function in ( |2.45| ). Both of these results relate the corresponding anomaly to 
a charge asymmetry. 



Appendix D. Energy one point functions in theories with a gravity dual 

In this appendix we present the calculation of energy one point functions for states 
created by current operators and the stress energy operator. 

D.l. One point function of the energy with a current source 

We wish to compute the contributions to the energy one point function ( |2.3U| ) for a 
state created by a current operator at strong coupling. The AdS/CFT dictionary says 
we need to compute the bulk three point function between two bulk photons and the 
graviton. The bulk action ( |3.18| ) contains two terms. The first term contributes also to the 
current two point function while the second term in ( |3.18| ) does not. Thus, the first term 

65 



contributes to the part of the energy one point function which is determined by the Ward 
identity and the second one to the angular dependent term which is not fixed by the Ward 
identity. In principle, the first term could also contribute to the angularly dependent part, 
but we have argued, based on the results for Af = 4 SYM that it contributes only to the 



constant part. Thus in order to compute the coefficient 02 in (|2.30| ) we need to compute 



the ratio of the contribution of the second term in (|3.18|) and the contribution of the first 
term in ( |3.18| ). 

Since the graviton is localized in W + and the photon is localized in the other transverse 
directions if the state has definite four dimensional momentum, we can approximate the 
computation as a flat space computation. In this particular case, this approximation 
will be exact, but that will not be the case when we discuss the three graviton vertex. 
Thus, we evaluate the vertex expanding the flat space action, but we will insert the AdS 
wavefunctions for the external states. 

In flat space our coordinates are (x + , x~ , x 1,2,3 ). The metric is ds 2 = —dx + dx~ + 
dx l dx l (latin indices % and j go from 1 to 3). 

We want to collect terms that are of first order in the perturbation h. There are 
two such terms, one per factor of g^ v in the action. The determinant g does not receive 
corrections and yfg = Our perturbation is of the form h = h+- ^{x l : n l )5(x + )(dx + ) 2 . 
We made explicit the fact that h ++ depends on the transverse coordinates and on a unit 
vector n % in the transverse space that represents the position of our calorimeter. 

Therefore, we want to calculate 

Si = ~ J dX+dX ~ d3X 2h ++ F+ i F^g ij = -^L J dx + dx-d 3 xh ++ F + *F+3g tJ (D.l) 

Notice that the contraction of Fs is restricted to the 3 dimensional transverse space as the 
metric element g-- is zero. We can do the x + integral easily as h++ is localized in this 
direction. Also, we will use the fact that the wave function of the photon is localized in 
the transverse space ( |3.15|) . We represent this fact by writing 

F +i F +j gij = a(x+, x-)S 3 (x l )e l e j giJ = a(x+, x~)5 3 {x l ) (D.2) 

where e l represents the (normalized) polarization of the photon. Notice that we choose the 
polarization in the transverse directions. Therefore F +l ~ d + A l . Using these facts and 
performing the integrals we get 

Si = — 3_ / dx+dx-d 3 xh ++ F +l F +J gil = — ^/i ++ (0,n l ) / cfaTa(0,ar) (D.3) 
4g 2 J 4g 2 J 
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In fact, h++ evaluated at x l = does not depend on n\ see ( |3.10| ) at W l = 0. Therefore, 
this term does not contribute to the angular dependence of the correlation function. 

Let us now look at the other term in ( |3.18j ). We first need to compute the Weyl 
tensor. Let us start with the Riemann tensor. This tensor has terms that go as \d 2 g and 
terms that go as gTT. Since we are in a flat space background only the first type of term 
contributes. This yields 

R +i+j = ^d i d j h++ (DA) 

All other terms are given by symmetry properties ( i.e. R+i+j = —Ri++j = —R + ij + = 
Ri + j + ) or vanish. The Weyl tensor also contain terms of the form \g\ u RnK- But, 

R^ = g Xp Rx^ — R++ = g ij R+i+j = \g ii d i d i h++ (D.5) 

We see there is only one non vanishing term that is proportional to the laplacian inside 
the transverse space. There are also terms proportional to the Ricci scalar inside the Weyl 
tensor, but we can see that these vanish in our case. The Weyl tensor is, then, given by 

C+ i+j = ~ (didj - \gijd k d^j h++ (D.6) 

The other components either vanish or are given by symmetry properties. There are four 
possible positions for the two plus signs, so we will have four terms in the second term in 
( |3.18|) (we are also using symmetry properties of F +l ). That is 



c a f dx+dx d3x Ar r+J - 

a 



Once again we can perform the integrals to obtain 



J dx + dx-d 3 xF +i F +J [d^j - ^gijd k d k ^j h++ 



(D.7) 



S2 = Udj - l -g l3 d k d^j h++(x\rf)\ x ^ei J dx-a(0,x~) (Di 



We are interested in the quotient between the angularly dependent term ( p.£| ) and the 
spherically symmetric term ( p.3| ) 



4a (didj - \gijd k d k ) h++(x\' 



\x i =Q 



M| h++(0,n*) 
We use the explicit form of the perturbation (|3.10| ) and we get the result 



(D.9) 



4 dS = -48^ (D.10) 



This gives the gravity result for the anisotropic part of the one point function (|2.30|) of a 
state produced by a current. 
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D.2. One point function of the energy with a stress tensor source 

Now we want to repeat this calculation for the case where we have the stress tensor as 
a source. In this case we need to consider 3 graviton interactions. There are 3 operators 
that contribute to this vertex. A natural parametrization is given by the action ( |3.24| ). 

We will do first the calculation in a flat space background. As in the computation 
we did above we will get derivatives acting on the perturbation h. When we go to the 
AdS background we could get terms involving the background curvature. Such terms are 
isotropic and will not contribute to the terms that have maximal angular momentum. But 
they do give contributions to the terms that have smaller values of the angular momentum. 
The computations we do here give only the leading contribution for ti and £4 in ( p.37|) . 



We start from the action in (|3.24|) and we expand each term to cubic order. We focus on 
terms with highest angular momentum in the transverse dimensions (see We use the 

fact that we need one of the metric perturbations to be h + + while the other two only have 
purely transverse indices. We find 

R = - ^ ++ ^ 1) (a+) 2 ^ 2) , R^aR^ = -28^+^(8+)^ 
R tlvS<T R Safn R f f v = -6d l d J d k d e h ++ h^ 1) (d + ) 2 h k ^ ) 

Notice that expanding the determinant of the metric in the action does not contribute to 
the three point function. If we now use that the wave function is going to be of the form 

fr(i)# +2 ^) = P( x+ i x-)5 3 (xy j e M (D.12) 

we can perform the integrals and calculate the quotients of the contribution to the three 
point function. After taking the derivatives and evaluating at x = we obtain the ratios 

^ 48 W (D ' 13) 
^ ii20 W (D ' 14) 

Due to the issues we discussed above, there are terms contributing to £2 which are of first 
order in 72/ (RAdsMpi) 41 , coming from the term with six derivatives in the action, which 
we neglected compared to the contribution of the four derivative terms. These formulas 
are also valid only to first other in the 7$. 
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